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A. BACKGROUND AND HISTORY 

Sergei Gorshkov, Admiral of the Soviet Fleet, has 
@iecenuvyea that the end of World War Two marked, "... 
the start of a military-technical revolution which in scope 
pm@emaepth transcended all the reforms and transmissions 
which had previously occurred in the armies and fleets of 
the world." [Ref. 1] 

Much has been made of the tactical and geo-political 
mepeeeatcions Of that revolution. Less popularly publicized, 
but equally radical are the changes that the developments in 
technology have spawned in the processes of ship design and 
memset ruction. 

iireOouehout most Of history, war at sea has been a 
process of trying to set one's opponent afire, board him, or 
ieee hunkS of Stone or metal at him across relatively short 
distances. All of these methods had the distinguishing 
Characteristic of being essentially disabling tactics, 
unless a ship was burned or perforated to the point of 
Sinking (and they often were not), damage was generally 
largely superficial and the ship could be refitted to eee 
service again. Witness the many vessels that have been 


mmeomt On botnm Sides of a contlict and the tradition of 


dt 


expecting naval personnel to supplement their pay with prize 
money. 

Two events occurred at or around the beginning of the 
twentieth century to change all of that: the invention of 
high explosives and the advent of the submarine. 

Initially, the submarine was primarily a commerce 
Peace Or It was too small CO Carry 4& prize Crew anmemenoe 
lightly armed to effectively combat a conventional warship. 
During World War One, most battles between submarines and 
conventional ships were fought on the surface. Submarines 
would use their stealth to sneak-up upon unsuspecting mer- 
chantmen, surface and then attempt to sink them with their 
deck cannons. Torpedoes were largely unreliable and the 
Submarines themselves were at risk if fast armed vessels 
were present. Convoying the merchants with Seis agile 
men-of-war was, briefly, an effective means of countering 
the submarine threat. However, as torpedoes became more 
effective so, too, did the submarine. For the convoy escort 
vessels it was no longer sufficient to have to wait for a 
Submarine to surface in order wo atiack ea bo effective 
Standoff weapon was needed. 

That standoff weapon was the depth charge. First intro- 
duced in World War One, the depth charge came of age in the 
second World War; and with it a new era in ship design was 
ushered into existence. For the principle of this new 
weapon was not one Of punching a hole in 42° Ship seagil oae 
then exploding internally; rather, high explosives permittee 
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generation of intense pressure waves in the sea that could 
break a vessel's back or rupture hull fittings whether the 
Vessel's skin Was punctured or not. Simply increasing the 
thickness of plating in discrete areas of vulnerability was 
no longer sufficient. Design for strength and resiliency 
under unequal pressure loadings was required. Moreover, 
ifmies requirement applied to attacker as well as attacked. 
Although the attacking surface vessel had the advantage of 
weet C LO ruil irom the immediate scene of the explosion 
and the effects were not compounded by as great a pre- 
existing static head of water, these weapons were not direc- 
tional and could be expected to damage all within their 
sphere. 

The issue of a vessel's ability to withstand an under- 
water explosion, be it submarine or surface ship, became 
more pressing as the explosive devices became increasingly 
powerful and sophisticated. Not only did depth charges 
become lethal over a wider area, but the advancing tech- 
Memoley Permitted mines, too, to grow from essentially 
contact weapons to devices that could endanger all types of 
Shipping from considerable depths. The issue has become 
particularly crucial with the development of nuclear weapons 
which are so powerful that no device need contact a ship or 
eeoup Of Ships to cause rampant destruction. 

The response of the U.S. Navy was to require that ships 


ama their Critical components be designed and built to 


iS 


withstand specified levels of shock loading. these re= 


quirements included (and continue to include) actual shoexk | 
testing of the first ship of a class and many crucial Sjeee7- 
on the component level. 

During the immediate post-World War Two period, many of 
these requirements were largely tentative and experimental 
By the early 1960's, however, sufficient experience had been 
gained to allow design requirements to be set forth | 
(Ref. 2]. These specifications established specific guide- 
lines for contractors to meet based upon type of equipment, 
tyve of vessel in which the equipment was to be used, and 
location and method of mounting within the vessel. AS with 
the ships themselves, actual shock testing of the item was 
the primary means of establishing compliance; aS 1t remains 
© Glave 

Performing a shock test upon production end predquer. 
particularly when those items are as large as an aircraft 
Carrier Can have Some ObV1IOUS Geavoj2e<- - it 1s costly waa 
1S dangerous, and the resales Can gee ibaa 
Moreover, there can be longterm structural damage which 
1S neither immediately apparent nor correctable, but which 
does not by itself constitute non-compliance of contract 
terms. Moreover, some items simply do not lend themselves 


T® ShGCK Test aaa. 
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What has been needed is a method of analyzing the design 
itself in order to achieve a numerical aoproximation of the 
end product's response to shock. 

This need has coalesced into the development of two 
computer codes; EPSA (Elasto Plastic Shell Analysis), 
produced by Weidlinger Associates, and USA-STAGS (Underwater 
ereck AnalysiS—--structural Analysis of General Shells) 


produced by Lockheed Missiles and Spacecraft, Inc. 


Eee OSI RCTIVES 

This study is the second in an ongoing series sponsored 
by the Defense Nuclear Agency (DNA) into the effects of 
underwater explosions upon ship hulls. The intent of these 
Studies has been to test the applicability of the two previ- 
ously Moneioned Structural response, underwater shock anal- 
YVS1S computer codes to naval design and analysis requirements. 
Under the umbrella of general avpnlicability, focus has 


been volaced upon two defined goals: 


1. Use of the codes to gain insights into large deflec- 
tion elastic-plastic responses of a Submerged struc- 
Pee moe jeclegd bOubLransilent acoustic shock loading; 
with especial emphasis upon stiffener tripping 


ohenomena. 

2. Performance of underwater explosion testing of the 
Structure to validate the results predicted by the 
code. 


To achieve these goals, flat plate models geometrically 
Similar to a ship's stiffened hull structure have been 


Studied. These plates were air backed to simulate general 


Lo 


conditions that could be Expecteq anaes plating; ata 
were subjected to shock loadings that produced deformations 
well into the plastic regime in Order to ensure Stitiewes 
tripoing. Existing information concerning underwater Crome 
Sion phenomena was employed in an attempt to reduce experi- 
mental uncertainties. 

The first study in this series was that conducted by am 
Thomas Rentz, USN in fulfillment of thesis requirements at 
the Naval Postgraduate School |Ref. 3] and |Ret yee 
that study he developed the following experimental strategy 
to study the EPSA code: 

An attempt to validate the code was conducted in two 
phases: Using pressure approximations based upon empiri- 
cally determined formulae, the code was first used as a 
pre-test predictor of test results. Once the underwater 
explosion testing had been conducted, pressure histories 
derived from that testing were inbut into the code we. 
plate strains predicted by the code in both its presShot me 
post-shot capacities were then compared te "26 pus mao 
histories from testing. Plate deformations and boundary 
Stresses predicted by the code were compared to observed 
physical responses. 

In this study, the code to be tested is USA-STAGS. aa 
an effort to maintain continuity with the vwearlicr study aaa 
the series, the same experimental strategy will be followed. 
The test geometry and flat plate model will also be 
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Metained. Conclusions ees acl ae aie —awriaey Ot Uan—-o AGS “tO 
predict the resvonses of a simple model to underwater explo- 
Sions will be drawn; and, perhaps of even greater impor- 
meme], Cuidelines will be established for uSing the 
USA-STAGS code and conducting underwater explosion testing 


Memmriucure studies. 


nee 


Il. UNDERWATER Bx Pi@s Tel eee on. 


Ase THE PRESSURES Oyo r 

Any conversion of matter, whether chemical or nuclear, 
which results in the very rapid pEoducme of large quant 
ties of gas at very high temperature can be categorized 
within the broad species “explosion.” Jlypically > {¢6. = 
tary high explosives such as TNT the pressure and tempera- 
ture within the gases produced by an exnlosion are on the 
order of 50,000 atmospheres and 3000 degrees Celsius 
[Ref. 5: p. 3]. For nuclear explosions the initial tempemie 
ture 1s on the order of a million degrees Celsius and the 
pressure, which is dependent upon the yield 1S similarly 
orders of magnitude Diener Geom tape on a conventional 
explosive [Ref. 6]. 

Clearly, temperatures and pressures of the levels 
described must dissipate high levels of energy through any 
Surrounding medium. In the case of water, the initial mani- 
festation of this disturbance is an intense compression 
created pressure wave propagating radially outward from the 
charge; followed nearly instantaneously by a displacement of 
ig) Ll In the immediate vicinity of the charge, the 
velocity of propagation of the pressure wave 1s Severa 
times the speed of sound (about 5000 feet per second in 
water) but approaches the value almost immediately and 
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therefore may be apprehended as an essentially acoustic 
Phenomenon. As the water begins to be displaced, the ga 
bubble created by the explosion expands; therebdy decreasing 
its internal pressure and correspondingly the pressure on 
the constraining fluid. Viewed on a pressure-time curve the 
imental rise in pressure 1S so nearly vertical as to appear 


Peeeoomncvinuous, while the pressure relief is characterized by 


mamapproximately exponential decay (Fig. 2.1) [Ref. 3: pp. 
4-7]. 
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MmeemG@e 2,1 simple Exponential Approximation of Incident 
Shock Wave 


The actual magnitudes of the pressure pulses and the rate 
of exponential decay have been found through experiment to be 


tS 


functions of the type of explosive, the charge weight, and 
the standoff distance.* The equations for the pressure and 
decay constant were originally determined by Robert H. 
Cole and were presented in his work, Underwater Explosions 
[Ref. 5]. These have Since been Simplified to use a series 
of empirically derived coefficients (@evelope4 b= 
Price: | Rei. sie 

The pressure profile can be expressed by the equation 


t, - t 


P(t) =P exp (eqn 2.1) 


where ae is the initial (greatest) pressure of the shock 

wave, 8@ 1s the decay constant describing the exponentias 

decay, andt - ts is the time elapsed from the arriva lee 
tie shock im omer 


The initial pressure, Poo and the decay constanim 


8, are dependent largely upon the type of explosive and 





the weight of the charge. These may be expressed 
Ie es 
W e 
= SS 2 
La cs ( a (eqn 252m 
and, 
gS NS 
_ «ge! 9(E ; z (eqn 2.3) 


ee standoff distance 1S deiined as tae distance fren 
the charge to the nearest point on the target. 


ZA 


~ 


where a Peooumescweer Saulare Inch and 3 in milli- 
seconds. or Ko» Ad VS are tne explosive dependent values 
(Table I). W and R are the charge weight (in pounds) and 


Standoff distance (in feet), respectively. 


Table I 


Exemesiveweonstants I 


Explosive coat ney s Sab PENTOLITE 
KS Pez IOS 22347 .6 24589 
Ay tees 1.144 dlipeael gS) 4 
Ko sai) oo 20S 2 
Ay leo -.247 et 


Eee SURFACE CUTOFF AND BOTTOM REFLECTION 

Shock waves vropagate radially from the charge. Thus, 
for a point source in an infinite medium, the wave could be 
expected to travel as an ever expanding sphere until damped 
by the fluid. For small to moderately sized charges deto- 
nated well below the surface in great depths of water, this 
Mewar Satistactory analogy. For explosions nearer the 
Surface or in shoal water, however, there are secondary 
effects for which there must be an accounting. 

When a shock wave reaches the surface of the water, the 
PemOmine air is of Sufficiently lesser density to be inca- 
pable of supporting the shock wave. The result is that the 


effective pressure at the interface must be zero in 


a1 


compliance with the laws of continuity. This implies st Gae 
a negative shock wave of magnitude equal to the original 
incident wave must be propagated downward through the 
medium to satisfy equilibrium conditions. This negative 
wave will travel through the water at an angle to the 
Surface that is equal to that Gf the incident Wave ween 

as the negative reflection image from a mirror continues at 
an angle equal to the angle of incidence. 

Remembering that these waves are radiating spherically 
Outward from the point of inception, it is Clear gna yee 
the incident wave and the generated surface wave will even- 
tually reach the ship, submarine, or other target object. 

It is likewise clear that, with its longer path, the surface 
wave must reach the target at some time later than the inci- 
dent; which travels only the shortest, Or Standoii, diocese 
between charge and target. The net effect is additive. 

When the surface wave reaches the target, the initial shock 
front will have passed and the incident wave will be at some 
point in its exponential decay. The negative surface wave 
will then "chop off" the tail Of themmnecidenue 2 [tf tie 
spatial geometry of the shot is such that the magnitude of 
the surface wave is greater than the remaining magnitude of 
the incident wave, a region of negative pressure will be 
Creaved Ct 1s 32a If this negative is less than 

the pressure required to keep the fluid in a liquid State, 
the fluid will "flash” inte Vapor ey tte Ne tO et ete eee 


Ze 





Beteryred tO as surface cutoff, and the negative pressure as 


mk CAaVitation. 
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Figure 2.2 Exponential Approximation with Surface Cutoff 


Pi tNeer eiiect Of bulk cavitation is, Sac Se lie not 
Meeerculariy harmful. The structure simply experiences 
Essentially no pressure while it is in the cavitated ie ce Ouae: 
When, however, the static head of water above the region 
Se -omes the Cavitation forces, the region closes Suddenly 
Bemecenerate a destructive reloading upon the structure. 

To further complicate matters, if the shot occurs in 
Shoal waters or at Gepths mear the bottom a second © Wake enie 7 


is also present. Unlike the surface, in this case the 


2 


incident wave Ws eines aS ne eneee medium of greater density 
than the fluid. The reflected wave will therefore be 
positive. The actual magnitude of this wave depends upon 
the physical characteristics of the bottom mater. ae 
effects of this bottom reflected wave are also additive 
and will generally reach the target at some time after 
the initial shock front has passed; altUhouegn thew aim 
time of arrival of this secondary front will one p22 
be dependent upon the bottom conditions and the speed at 
which the wave travels through the bottom medium. Addi- 
tionally, in shoal waters, that same bottom reflected wave 
will reach the surface; creating itS OWN Suni ace ecm 

The net effect upon the target 1S thus One Of ane 
dent shock wave; which, if the charge and target are suffi- 
ciently near the eee will be cutoff by the negative 
Surface wave. This, in turn, is, [ol lowed ey (howe e pec 
reflected wave which also experiences an exponential decay, 
and which may also be cutoff by its own generated surface 
Wey Camry Zee = 

Analytically, the easiest way to model the eGitecrts ues 
incident, surface, and bottom waves upon the target is as 
has been represented in Figure 2.4, with three separate 
images. 

The incident and surface images are relatively easy to 
depict. The incident image, May of course represents the 
charge itself and occupies the same location in the spatial 
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Figure 2.4 Three Image Presentation of Shock Waves 
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seometry. The Surface image, Me 1S or the same magnitude 
as the incident image but of an onposite sign. It lies on a 
Mine drawn vertically through the charge at a distance above 
the surface equal to the depth of the charge below. 


The surface cutoff pressure may therefore be ‘computed 
(eqn 2.4) 


where Re is now the distance from the surface image to 
the target. The time of cutoff after the passage of the 
initial shock front is thus the time difference for the wave 


to cover the greater distance, or 
- (eqn 2.5) 


where c is the speed of sound in the fluid. The cutoff 
pressure will then decay in the Same exponential manner as 
the incident pressure. . 

As may be conceived from the preceding discussion of the 
Boeewom reitlection; the location of the bottom image is not 
Semovcraightforward. It may be viewed as initiating along 
the Same vertical line as do the incident and surface 
images, and at a depth below the fluid-bottom interface 
equal to the distance of the charge above. Here, however, 
the similarities to the surface image end. The bottom image 
cannot be considered a point charge because of the tendency 
of the pressure wave to "smear" across the bottom. The 
magnitude and speed of the wave, too, are dependent upon 


ET 


bottom characteristics. A mudady DGtTOmeior 10S Gee eee 
tend to reflect a diffuse, low magnitude wave travelling at 
or about the sveed of sound in water. A rock bottom on the 
Other hand will produce sharp, high magnitude pulses waver 
will tend to radiate over a distance due ‘to the higher speed 
of travel of the shock wave in the bottom medium than ine 
fluid. Due to this dependence Upon DoOllom Conair 
determination of bottom reflection effects do not lend then 
selves to general predictive studies and are often disre- 


garded, although’ they can beget on it ieanee 


GC: BUBBLE BEE rer 

The foregoing discussion of the development of the pres- 
Sure wave and its impact upon the target assumes that a 
Single, clean wavefront is produced by the explosion. In 
point of fact, such 1s rarely the case. When the charge 
explodes, it fLirst creates, a0 ren intense Shock» £renme 
which 1S the initial incident shock wave. It also vrodug 
a bubble of very hot gasses at very high pressure. 

The bubble expands until 1t reaches pressure equi 1 ioe 
with the surrounding water. When it does so, a portion of 
the energy within the bubble is released as a new shock 
front. This release of energy then permits the momentum of 
the water displaced to assert itself and the bubble is 
collapsed, compressing the gases within. The displaced 


water can be conceived as attempting to regain its original 


ao 





pre-bubble condition. As the water rushes in, it compresses 
the gases beyond the equilibrium condition and the bubble 
again expands, again releasing a portion of its energy as a 
shock wave, albeit a smaller one than the preceding (Fig. 
2.9). This process has been known to produce aS many as ten 
Smeenmiticant pulses [Ref. 53: pp. 8, 9], which, though weaker 
than the initial shock front, can still have a pronounced 
effect upon a structure which has already been loaded by the 


earlier pulses. 
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Figure 2.5 Multiple Expansions of the dot Gas Bubble 
feet . of Dp. 7) 
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The time and radius of the bubble at the first pulse 


has, again, been determined empirically, the maximum radius 


as 
wil 3s 
Amax Dean (eqn 2ige 
at time, 
wi/ 3 
cia ‘(a5 (eqn 2.7) 
where ae is expressed in feet and t in seconds. Wis 


Once more the weight of the charge, Ke and Ky are explosive 
devendent constants (Table II), and D is the depth of the 


charge in feet. 


Table I] 


EXp LOS1Ve P@@ne can ts ai 


Explosive ABX-1 TNE PENTOL Ti 
Ke 14.14 iP etoea 12 Oe 
Ky Aa tous 4.268 4. 3a 


A second and distinct phenomenon is bubble migration; so 
called because of the tendency of the bubble to migrate 
toward any nearby mass. When the bubble contacts the mass, 
it collapses and begins a series of rapid pulsations. The 
net effect is Similar to that of the already described 
bubble pulse; a series of reloadings upon the structure. 


The rapidity of the pulsations of the collapsed bubble and 


30 





mietr increased localization tend, however, to magnify their 


severity. 


See PLICATIONS OF THEORY TO EXPERIMENT DESIGN 

Developing a basic understanding of the primary 
phenomena involved in underwater explosion theory is 
crucial to developing a successful underwater test program 
for many reasons. Two are of particular interest here, 
however: 

mest, ana most Obviously, equations 2.1 through 2.3 
allow initial estimates of the peak pressure and pressure 
history to which a test structure 1S to be exposed. These 
are then input to the computer code so that preshot approxi- 
mations of the structural responses may be made; thereby 
emeacing guidelines for calibrating gages, determining 
charge weights, etc. 

PecoOna, ene design Of the test itself is a direct conse- 
quence of the phenomena being observed. As is readily 
apparent from the earlier sections of this chapter, the 
loadings imposed upon a structure by an underwater explosion 
can be extremely complex. The incident wave may be a 
Simple exponential decay dependent only upon charge and 
mmeamaott, but the many secondary effects are subject to 
variables which cannot be readily controlled experimentally. 

The contributions of the secondary effects are signifi- 
cant and must eventually be studied. However, during these 
early phases of investigation into underwater explosion 


oe 


effects a determined effort has been made to simplify the 
problem as much as possible. 
It has therefore been required that: 


1. The model to be studied be aS eae 15 cone 
tent with stated objectives. 


2. otudy be limited only to the ancident shock “wave 

These two requirements are interrelated and are tied 
intrinsically £o the nature Of underwater exp learn 
phenomena. AS has been seen earlier, the incident shock may 
be approximated as a plane wave emanating from the charge 
and striking the panel. In the absence of secondary 
effects, the responses being studied here may thus be 
limited to those of the test panel alone. If, On the Otitem 
hand, secondary effects are included, more than a single 
point source is involved and the interactions of any SUDDpGia 
or backing structures with the fluid medium must 2c 
considered. 

Fortunately, for the purposes of this study, suviaa. 
cutoff and bottom reflection are later time responses which 
can be neglected if study is confined to @he fimst tew 
milliseconds after detonation. Bubslemiieees en 
effects can be reduced greatly or eliminated through careful 
attention to the test fpeometnry. 

If.an exploSlOneoccurs ssi fase © ae ee ee 
that the initial bubble radius is greater than the charge 


depth, the hot gases forming the bubble will vent into the 
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atmosphere and be dissipated. This venting iS instantaneous 
and, if occurring at the proper depths, prevents any further 
beroole Dulsation or migration. 

Bulk cavitation, too, is dependent upon the test geom- 
meee Lhe region Of Cavitation tends to expand horizontally 
away from the line drawn between source and surface images 
meee 2.0). tte actual Shape and dimensions of the cavita- 


“ 


ifemereeion Will vary with charge size, type, and depth, but 
The area at depths below the charge can generally be 


expected to remain free of cavitation effects. 





IMAGE Nie 
SOURCE 77 ~~ 








CLOSURE yy ty, 
. Lip“ yyy *® “UPPER CaviTATiON 
= Ly NY wo ENVELOPE 
See My 
=. 
Wie oO “ =o 
SOURCE XS. Si ase aa LOWER CAVITATION 


ENVELOPE 
BHOCKX 


| 

| 

| 

< 
> 
| 





PiwC mer omeneofton Or Bulk Cavitation 
(From Weidlinger Associates CUE Code. ) 


oo 


Bulk cavitation and bubble effects Gan (thereto w ae 
reduced Sigenit team! 3 eeu Sane entirely it the test 
peometry is such that the charge is located directly above 
the plate at a depth to cause bubble venting. 

In designing the test platform [Ref. 3: pp. 19-28] a 
series of iterations were run uSing the EPSA code to deter- 
mine a charge weight and standoff distance that would 
produce maximum test panel deflections of approximately four 
plate thickness. - From these runs it was determined that 
the ideal charge weight and standoff were 8 pounds of TNT at 
Set Ce u., 

The maximum bubble radius is expressed as a function of 
charge weight and depth (Eqn. 2.6). Since it has been Eom 
that bubble venting occurs when the ratio of charge depthmim. 
maximum radius is less than ase choosing a charge depia 
is a fairly straightforward procedure once a determination 


of the desired ratio has been made. 


See value was recommended by Weidlinger Associates as 


One which experience indicated would ensure stiffener 
Ue aoe. 


2 ico ratio below 0.75 will cause venting. Venting by 
its very nature will duct a portion of the explosion energy 
into the atmosphere, thereby decreasing the maximum pressure 
of the shock wave. As a rule of thumb, this decrease in 
usually considered negligible until the ratio drops below 
6) air 
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Pesce locrmme a2 Ceptn tO radius ratio of 0.30 it was 
Pel, desrermined =emat a charge depth of four feet was 


meeropriate. Applying this to Fan. 2.6 yields 


gl/3 
A = (12.67), ——--——— 
max oe 33)1/8 
= fo Fea 
a iG 
lle ») 
Colt 0 


which assures venting. 

From the foregoing it can be seen that a test Senet cura= 
tion that will have the desired consequences of reducing 
bubble and bulk cavitation effects can be designed. For an 
eeeeme DOUnd TNT charge, this platform will have the charge 
centered vertically at a standoff distance of nine feet 
above the plate. The charge itself will be located at a 
Seem Of four teet below that surface. The actual design of 


the test platform will be discussed more fully in Chapter V. 


te 


Ill. THE FLAT PLATE MODEL 


As stated in the introduction, the OdjJective Gi time 
study has been to test the applicability of a computer code, 
USA-STAGS, to predicting the effects of an underwater explo- 
Sion upon ships! hulls. To do this, it 1s clearly neces ama 
that some basis for modelling simple Sections Of a type 
ship's grillage* be established. 

The scaling laws used here and in the earlier work of 
this series [Ref. 3] were those developed by Dr. Raymond P. 
Daddazio of Weidlinger Associates, Inc. 

In his work [Ref. 8], Dr. Daddazio has established tig 


dimensionless parameters, 8 and A, which ne Tre{eeS elon 


the "plate slenderness ratio," and the "longitudinal stif- 
fener slenderness ratio,'' respectively. These are defined 
0 
om cae (eqn Gaim 
it E 
and 
a oy 
CS See 2 
r ee 3 Ceqn 3.2) 


=e ship's hull Or supexmstrucitural plating (wey as occu 
ated stiffeners. 


SO 


—— Mategmereyiecld Sstrengi na 

fee > 1OUne S modulus 

a = spacing between transverse frames 

b = spacing between longitudinal stiffeners 

tC = plate thickness 

feee- radius of gyration of the longitudinal stiffener 


acting with an assumed effective width of plating 


Oo. where, 


i 


g¢ 


oa, 


WD 
mip 


ireoewol Samobing a variety Of Ship types, Dr. Daddazio 


Mumma that typical values for ships’ grillage fall into the 


range: 


vnc. 


Mie intent of this 


Semiea more cte lates has peen to use one 


Dasic model but with a variety of stiffener types so that 


the responses of different stiffener geometries could be 


Sealed. 


Mimo delmestoOl Iehned f@mahet. 3 was a 3/716 inch thick 


test panel 18 inches long by 12 inches wide milled into the 


Peer Ole Evo ineh thick 6GOGl1-T6 aluminum blank measuring 


a7 


33 inches Dy 27 Mimemes. Two stiffeners with rectangular 
cross sections, 3/16 inch thick by one inch deep were mounted 
transversely across the test panel in Such a manner as to be 
symmetric to both the transverse and longitudinal axes (Fig. 
3.1). The entire assembly was machined as a single unit to 
avoid the inconsistencies of welding. A radius was cut at 
all corners and the whole polished to avoid stress concen- 
Craton POE ee 

To apply Equation 3.1 and Equation 5.2 Ghee eaaeee oe 
stiffeners were considered to be longitudinal stiffening 
members and the longer boundaries of the panel to be 
transverse frame members. Thus, the dimensions of the plate 
are, a = 12 inches, b = 6 inches, and t = 3/16 inches 7) Ue 
the standard reference values for 6061-T6 aluminum, it was 


found that, 16 90 ats ecko Vy ae eon acne: (Ref. 3] 


B= 2.02 
ae! 
XN = 0.9 


which fall within the accenotable ranges. 

The initial plate used in this study was identical to 
that used in Ref. 3 with the exception that the stiffener 
cross section used was a tee rather than a rectangle (Fig. 
See 

Since the spacing of the stiffeners in this tee- 
stiffened plate was identical to that in the original plates 
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Pigure 3.1 Test Panel with Transverse Rectangular Stiffeners 
[_Remee@o: . pb. 13] 
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clearly the 3 values remained invariant. Given the 
changed geometry of the stiffener the effective radius of 
Pyration of the plate/stiffener assembly (k) could be 
expected to change, and with 1t the value of i. 

Before that value A could be determined, however, the 
@emensions of the stiffener itself first had to be 
eeeaolished. 

In an effort to ensure that stiffener responses were a 
meme t1On of stiffener design alone, the plastic moduli, oe 
of the two stiffener cross sections were set as close as 
possible to equal. For the original, rectangular-stiffened 
peace , 3 equaled 0.046875 in’. 

By a process of iteration, a tee-stiffener with flange 
dimensions 23/32 inch by 1/8 inch, and web dimensions 25/32 
maem OY 1/8 inch was fixed upon aS giving an elastic modulus, 
Z = 0.046170 ae Viiveeasor Veith Deane Omumirman tor the 
Meebangular stiffener (Fig. 3.3). 

Now, applying these dimensions to the grillage param- 
eter, A, we find that the neutral axis of the effective 


eeoitener (Fig. 3.4) is: 


7. i 
An 


Ore 22 OF in 
Mmeen moment Of inertia of the efifective stiffener is: 


= f 2 
log = ot + LAG. 


peOODscoean 
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=a Bees 2 Sais Neutral Axis: 





LALY, 
oa eo ae 
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32a 32 8) ee 16 
2S esl 23, 
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Oey 14.) sues 





Plastic Secrelon iodine. 


4 Z = Sum of the moments about the 
neutral axis of the cross 


? k— hb, sections above and below the 


| axis 
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feeetice 3.4 Eitfective Cross Section of Tee-Stiffend Plate 
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which gives a radius of curvature: 


Log 
A 





sale 
2 


2 
[0.092562 in*/1.0256 in“] 


O's 30042 ite 


II 


Thus, for the tee~stiffened plate, 


tlh 


Bee 40000 
-300419T 35. 196 
= 0.804144 


which is well within the parameter range required for 
modelling ships' grillage. 

The second requirement levied upon these test plates was 
the ability TOmdemonst rates St 1 ie geist ase ereni 

Stiffener tripping has been identified as the lateral- 
ii sano ent instability of the stiffener as it becomes 
Suddenly unstable and fails under load. It is identiziem by 
a characteristic warping and buckling of the stiffener. Tf 
1S Significant that the continued resistance of a stiffened 
structure to deformation under loading 1s diramartieaes. 
reduced after the stiffeners have tripped. In earlier 
studies of impulsive loadings upon ships’ grillage;, trippmee 
has been identified as ". . .a primary ductile failure mode 
for ship strucwure ss Rel oes. Oren 

Previous observations have indicated that tripping is 


most readily initiated when an axially compressive loading 
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memepolied to the top fibers of the stiffener. Recognizing 
that a flat plate will tend to deform away from the source 
Speen Explosion, the plates in this series of studies have 
been oriented with the stiffener sides outward into the 
ioe Cd medium.° 

AS stated earlier in this chapter, the original intent 
in this series of studies was to maintain an essentially 
fixed plate model, varying only the stiffener geometries. 
Unfortunately, the results of the first two explosion tests 
(on the original rectangular-stiffened plate [Ref. 3], and 
On the tee-stiffened plate here) were not satisfactory in 
view of the stated objectives. The ~wO Uraliover so mounted 
stiffeners were too short and the entire assembly too rigid 
memerttectivelyecemonstrare the desired stiffener trappine. 
Mie 1OCK wave essentua@bly ‘punched through the panel 
material, shearing it at the fixed boundaries and at the 
base of the aeons CE Gener). 

Mee PesurtUS Gi ethne two initial tests andicated that a 
@esien change was in order. The third plate of the series 


Meme second of this study) was a simplified model to permit 


5 ; : 
This is hardly the expected orientation on a surface 
vessel, but 1S common on many Submarine pressure hulls. 


Se eeenouen Meu Casy 16 See dm Proprers oy) noteall of the 
damage to the tee-stiffened plate was done by the incident 
Shock wave. Apparently, when the bottom reflection hit the 
pete and its backing structure, it was forced upward into 
its Own Supporting cables. Much of the damage to the stif- 
meV mweomli OAareclLcular seems to have been caused by the 
Herons Of that cable. 
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efoser investigation of the shock eifects themselves, 
without attempting either to retain a set plastic modulus 
Se tO approximate ships' grillage. 

Tne result of this simplification was a thickening of 
Ehe base panel to ¢ inch to increase resistance to shear. 
The stiffener was restored to the original rectangular 
@emrrouraction but the thickness was increased to = inch 
and the depth to 14 Ae VeS@me linen canr, i CHiLs 
attempt, a single, longitudinally oriented stiffener was 
used rather than the two transverse stiffeners in an effort 
to provide greater flexibility along the length of the 
pmeeatener (Fig. 3.6). 

Although the charge used to test the rectangular stif- 
fened plate was small (4 pound of TNT), and the correspond- 
ing plate deflections were far too slight for any firm 
conclusions to be drawn, it can be seen in Figures 3.7 and 
Peromthat tripping action appears to have begun to manifest 
Meee lr., ithe Single, longitudinally oriented stiffener seems 
meme Lne Superior model; and, while further testing is 
required at larger charge weights, its employment is tenta- 
tively recommended for future studies. 


° ‘This decision was made somewhat arbitrarily. Dimen- 


memes were Chosen, in part, to reduce the cost and simplify 
fren inineg the plate. it is interesting to note, however, 
miei with 6@ = 1.52 and 4 = 1.1, this plate is not far 
outside the acceptable grillage range. 
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fee GENERAL 

EeructuUral Analysis of General Shells (STAGS) 1s the 
finite element code first developed by Lockheed Missiles and 
Space Company, Inc. during the 1960's in response to an 
mreeial tasking by the Navy and the Air Force plus an inter- 
nally perceived need for a general structural analysis 
preogram. 

ime w"ourrent version Of STAGS is STAGS C-1, a dynamic 
analysis code with broad applications to the aerospace 
Mmeemstry which has also found application in ship design and 
engineering. A new version, STAGS 2000, iS now under develop- 
ment which will have an entirely changed data manager 
concept from that used in the C-l version. This change 
Should make the new code both easier to use and far more 
economical of the user's time. STAGS 2000 will appear 
Mmmemially aS a static analysis code, but will eventually 
exceed STAGS C-1's capabilities. 

STAGS C~-1 is supported by Computer Software Management 
and Information Center (COSMIC) under the sponsorship of the 
National Aeronautics and Space Administration (NASA). 

The USA-STAGS code iS an extension of the basic STAGS 
Pomnine sO Dermice Gyaluation of impulsive loadings upon 


Meo mesprUGCltures in a fluid medium. As such, Underwater 


oa 


Shock Analysis (USA) is a multimodute Toucine whiten 
inserted a an intermediary and adjunct to the program pre- 
analyzer (STAGS1) and analyzer (STAGS2). The STAGS?) Sloe 
combination can stand alone for some types of general struc- 
tural analysis. AS a result of itS add-on Nature, the Use 
code cannot be used in a stand alone mode, but it is not 
STAGS dependent. Versions have been created for use with 


Other finite element Structural analysis Soden 


Be THECRY 

The characteristics which make USA-STAGS unique lieu 
its applicability to a wide variety of problems and, Dai-taae 
ularly, in the individual techniques used to solve those 
problems. The broad, underlying theory behind those solu- 
tions 1S in itself common to most finite element codes. 

Briefly, the technique used by codes of this ore 1S "te 
solve the differential equations of motion fOr a4 Strucuma. 
in order to determine its responses to itS environment 
[Ref. 10: pp. 2-1 through 2-14]. By then ang tao 
constitutive equations (Hooke's law generalized to include 
anisotropic materials), the stresses in the governing equa- 
tions of motion may be expressed in terms of strains. These 
Strains may in turn be expressed as the first partials of 


the displacements with respect to the spatial coordinates 


As of this writing, USA has been linked to NASTRAN, 
GENSAM, SPAR, and DIAL. 


a2 


@ieouth application of the Kinematic relationships. an 
equivalency can be drawn between the motions of the struc- 
ture between any two instants in time ty and ty and the 


extrema of the energy functional sFacu - K)dt with respect 
moma! admissible variations of ee, energy U and the 
kinetic energy K which satisfy the geometric boundary condi- 
im~ems and COntLinuity conditions. The equilibrium equations 
Can then be derived from the first variation of the 

energy functional through partial integration with respect 
to time and the spatial variables. 

Miiemvarlatlonal approach provides a particularly effec- 
tive means for analyzing continuous systems because of the 
Meayeeit) Which it generates and accounts for natural boundary 
conditions [Ref. 11]. The disadvantage of the method is 
mere it 1S limited to conservative systems. However, by 
applying bie Deiiemere  OmmwrtlAlework woleaGco.i1S. able to 
greatly extend the classes of problems to which it can be 
soobi1ed. 

imaetne Second case; where inelastic deformations such as 
plastic strains or thermal expansions are encountered, STAGS 
ieeats them as “pseudo-loads." In this manner, by applying 
Meee pseude—-load cerms) tO the force side of the equation, 
an initial estimate of the inelastic deformation values can 
be made and, if convergence is obtained, the nonconservative 


problem can be solved by iteration upon a Series of 


conservative problems [Ref. 10: 002575) 2-3 |. ete ane 
of cases to which STAGS may now be a@oplied 1S Site ea. 
Figure 4.1. 

After being reduced to functions of displacement on ae 
the governing equations Of motion for a discretizea = 73. 


may be expressed in the form: 
. a ee E + a a \ a 
M_< 4+ D(x) 1 BC mega. Oo eames 


where x is the vector of displacement components, Me is 
the structural mass matrix, t is the vector of eGxternaa 
forces, and Ke is the (generally nonlinear) structural 
stiffness operator. The operators oe and Di include 
forces that are functions of structural deformation ey ma. 
deformation velocity, respectively. 

Equation 4.1 is the matrix equation thar yo! 4G. emee 
gos to find structural displacements. For [luivd7 seeuwee 
tural interactions the generalized force vector becomes a 
much more complex entity [{Ref. l2: pp. 2-1 throvwen 2-22 


For submerged structures excited by an acoust teswave 


the exciting forces are given by: 


f = -GA,(P, + Pg) (sqmaa4 25 


where Po and Po are the incident and scattered nodal pres- 


Sure vectors acting upon the Surface of the Structure at tie 
Cfluid=-sStructural interreaee= A 1S the diagonal area mat irae 
associated with elements in the fluid mesh, and G is the 


o4 
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matrix which relates the structural and fluid suriace 
forces. It can be seen from the development of tne pressure 
characteristics of an underwater Shock wave in Chaprter aa 
that for an underwater explosion the incident pressure 


vector Po is known. The scattered pressure vector Po on the 


other hand is dependent upon the configuration of the Sura 
tural surface and is therefore net enly not amediace., 
known but is also changing as the Strweture deme 

It is to solve this computational dilemma and thereby 


solve the combined equation of motion: 


MX + D(X) ~ Box) + K(X) = “GA .(P. + Pp 


Po) (eqn 4.3) 


that the USA code was created. 


The USA code finds P, by utilizing the Doubly ASymptGaae 


S 
Aoproximation (DAA) and a staggered solution procedure. 


The DAA is an approximate relation which approaches 
exactness at both the high-frequency (early-time) and low- 
frequency (late-time) limits to the passage of the shock 


i Pema It may be writteus 


M.P 


f + ‘Oe 


. ees = OCM Vg (eqn 44a 


where Vo is the vector of scattered-wave fluid-particle 


velocities normal tevthe Samet) Siete ate aoe 
Structural interraeer M 1s the £luid mass matrix for time 


fluid mesh at the interface, and 9 and c are the density and 


sound velocity Of Mime iii eK cee rl ely. 


Do 








Ween Llimtlea ana Structure are excited by an acoustic 
Weve, the incident and scattered fluid particle velocities 
V 
Oe 


kinematic compatability relation: 


and Vg) may be related to the structural response by the 


Ge = v. + v Ceqn 4.9) 


Bervine for v and substituting the resultant into equation 


4.4 yields 


is oem. (Gs — (eqn 4.6) 


MpP ai 


f och 


Ps 


S ie 


EFquations 4.3 and 4.6 are now each expressed in terms of the 
Single unknown P. and a solution may be obtained by a Stag- 


~ 


Sered iteration process between the two equations for each 


> 


time step. 


tee ORGANI ZATION 
As has already been stated, USA-STAGS is a modular code. 


memeclOonally, it may be divided into five basic components. 


fee tne structural preprocessor 

2. The fluid mass preprocessor 

3. The augmented matrix preprocessor 
4. The time integration processor 

Se lhe data postprocessor 


Each of these components is run individually in sequence 


utilizing the outputs from the previous modules. 


ae 


1. The Structural Preprocessor 
The structural preprocessor io tie Unlt sae 
assembles the structural mass and stiffness matrices, 
creates the finite element mesh, and otherwise provides an 
encoded description of. the structure's internal aie 
ternal degrees of freedom [Ref. 12: pp. 3-3, 3-4]. 
For USA-STAGS, the structural preprocessor iTS )7a 
STAGSL portion Of the SlAGS cede" tr are 
Input for STAGS is made to each of ten ~iunev prema 
subdivisions which correspond to the types of input thames 
Into. cach Theses) eGalv isle ees 
- Summary and Control Parameters 
~ Computational Strategy Parameters 
- Data Tables 
- Geometry and Material 
= Discretizacion 
- Discrete Stiffeners* 
= BOuUndary Condit homes 
- Loads* 
= QUGPMeG CenurorL 
- Finite Element Units* 
(Subdirectories designated with * are not required with all 
models. ) 
Input to each SubdLYreCUGry.) 1S ieee ss amor 


series of data records eensisStame Cigrealeend eae 


oo 





numerical data fields. This will be discussed in greater 
detail in the following section. 
Zee lL hew ruta siass Preprocessor 
The fluid mass preprocessor (FLUMAS) constructs the 
moc Mass Matrix for a structure submerged in a fluid. 
iees include [Ref. 12: pp. 3-4, 3-5] the global coordi- 
Memes fOr the structural geometry (Cinput from the structural 
meeprocessor), plus 
~ Fluid wet-surface mesh 
- Element definitions 
meiaverial properties Of the fluid 
eee. ct 1lOngor free Surface 
~ Model symmetries 
3. The Augmented Matrix Preprocessor 
Pic Valet ce Maem Ore wEOeessOr CAUGMAT) receives 
Mme@emmoutput from the structural and fluid mass processors to 
create the matrices required for solution of the augmented 
equations 4.3 and 4.6. 
4. The Time Integration Processor 
The functional heart of the USA-STAGS system, the 
meme winterrator processor (TIMINT) links back to and supple- 
Memus the structural processor STAGSZ to implement the stag- 
gered solution technique; creating the solution data files 
Over the time range specified. Input includes the data 


passed from AUGMAT plus 


ae 


- Spatial geometry of charge and target 


- Initial incident pressure and pressure wave decay 
characterise ies 


- Time integration information 
- Restart data 
TIMINT includes an internal poStprecessinye cape o ween 
Generally, however, this function is left to the following 
module. 
Oo. The Data Postprocessor 
USA-STAGS Provides a selection of three possible 


postprocessors: POSTP, POSTPR, and STAPL. Their uses will 


be discussed in the following section. 


De Uc LNGl ire meen e 

As has already been stated, the STAGS code is extensive 
in its concept; attempting to cover ina single code as 
many approaches to problem solution as can be linked in ; 
a Single contextual framework. 

It 18S beyond the scope of this work Gomeco 7-4 a 
possible aspects of the code; indeed within the stated capa- 
bilities of STAGS this approach may seem peripheral, being 
limited as it is to those aspects clearly related to under- 
water shock analysis of flat plate models. 

STAGS is designed to be used for any one of six types of 
analysis [Refi. 13: p22 oar 

- Linear stress analysis 
- Geometrically nonlinear elastic stress analysis 
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[ene lAst@e stress analysis, f@eometrically linear or 
nonlinear 


eoiturecabtonm DucKling analysis With linear or geometri- 
cally nonlinear stress state (elastic) 


- Small vibration analysis based on linear or geometri- 
cally nonlinear stress state (elastic) 


= Transient response analysis, linear or geometrically 
mombinear, elastic or inelastic 


In this study, only two of those capabilities were used: 

Oe Small vibrations (in a stress free state), and (b) 
Geometrically nonlinear transient response analysis. Of 
these, the second alone is truly applicable to USA-STAGS. 
The small vibration analysis was used as a check on the code 
Mlcemimplt Gata. Since the actual inputs for the two anal- 
yses vary but slightly, both shall be addressed together. 

Before that can be done, however, the codes must somehow 
be brought into some semblance of conformity; for they are 
Mememormally provided as a single unit. 

Perl COdeS are provided in tape form: STAGS from 
MeerecC, and USA directly from Lockheed. These are designed 
to be loaded onto the VAX 11/780 computer using the standard 
MOUNT MF: command.” 

muitos SUSFeStea That, at least during the early stages 
Pmeprocessing, each tape be assigned to a separate direc- 


mery, £Or Many OL the source files will be culled as the 


“See Appendix A. 


oak 


executable files Ca are created. The source files 
Supplied with each code are specific to the code until the 
final linkage, and their sheer number Can Cause cont US i ¢@iaee 
not approached methodically. 

The STAGS tape as provided contains all the necessary 
ingredients for the functioning codes except the executables 
themselves. These are left to be linked on the individual 
computer. The link commands for creation of these execu. 
bles are provided with the tape on the INFO.VAX file, which 
also contains a full listing of the tape's Contents aa 
commands that are of immediate interest are LSTAGS1.COM, 
LSTAGS2 .COM, LPOSTP.COM, and LSTAPL.COM. These Shou tai 
copied into the file and executed in the normal manner 
LOr Ta {progr ainsi. Gems eae Beene aie The result will be 
the executable Tlies Oru ne sae Spee maa mono 

Executables are provided with the USA tape; however, 


Since relinking may be required to adapt the code to local 


9 
filenames, dig k “commanGy Sales seme e vise pee cee The 


on the VAX, the asterisk represents a "wildcard" which 


may be substituted for any filename, filetype, or version 
number (the Gomplete 1istine for a 2116 oe ft bomen 
filetype, version number) whenever a general category is 
desired or the listing is to be repeated. 

This is done by simply typing "'@ command filename" and 
hitting the return key. 

Te USA tape also includes another particularly nie 
feature, the complete user's manual for each of the modules. 
This ensures that the mMOst current version 1S avai aoe 
and, more important, thewversicn of the Manta eee 
pertinent to that tape. A vital accompaniment since the 
Tapes are continually Sette wo dared. 
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command files to be used in constructing the USA executables 
meee CP eUeCONe  INKRALG. COM, LINKUSAS.COM, and LINKPOS.COM. 
Mmeeets the LINKUSAS.COM file which actually performs the 
marriage between the codes; allowing the integration/ 
Pameieotural analysis functions of STAGS to be coupled with 
the USA fluid processor. 

Once all of the executable files have been created and 
[Meoereod in a full run of the combined code; all files but 
phe executables may be deleted, thus realizing a considerable 
increase in computational Ce 

The basic strategy that will be followed here shall be 
to address USA-STAGS as a Single code broken into five 
different modules (see preceding section); each of which can 
be run independently in stepping stone fashion once neces- 
Sary input from earlier modules has been obtained. 

ii mi testemocwle Go be run is SIAGSI. The actual input 
for this includes all information needed for both STAGS1 and 
STAGS2 execution in a single entry file. 

ite first tilve daca records: A-1, B-1, B-2, B-3, and 
C-1 are common to all STAGS runs. The inputs as set forth 
Mmomerie COSMIC STAGS users guide [Ref. 13: pp. 3-1 through 
Peet) are largely self-explanatory but a few points should 
me noted. 


eo tal Space required for the executables alone is less 


mires) POlLOCKS, Compared tO the almost 13000 blocks for 
the executables plus source codes. 


ee 


First, IPOSTO, IPOST1L, and J[POST2 are applicable on ieee 
ties o LAGS post processors POSTP, and STAPL. If the Wserae 
interested only in the nodal displacements resulting from 
underwater shock these may be derived directly from either 
of the USA modules TIMINT or POSTPR. Thus IPOSTO, IPOST1, 
and IPOST2 may be set to zero thereby suppressing creation 
of model and solution data files (*.MOD and FOROZ2 Dail) 3 
realizing some savings in computer space. It is suggested 
that this be done in any case during carly model déevelopm aaa 
and while testing code inputs (as with the small vibrat ioe 
ieael Vs Sor If, however, strains are required; these are 
computed only within the STAGS postprocessors and I[POSTO 
must be set to 1. Considerable amounts of computer T inieue 
Space can be used in processing a large model; time wawea 
would be wasted if this is not done. Similarly, IPOST1 and 
IPOST2 should be set to a non-zero value if displacement and 
stress plots are desired at the IPOSTist and IPOST2nd load 
con 

second, it should also be noted well that IFLU musta 
Set to 1 if underwater shock analysis using USA is to be 
undertaken. This will cause the FOROO3.DAT and FOROO4.DAT 
files required for USA to be created. These files are not 

ae is required only when using STAPL. POSTP is 


capable of calculating displacements and stresses 
internally. 
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unduly lengthy, but for the sake of economy should be 
Suppressed for non-USA runs. 
The C-1 data record is meaningless for small vibration 


analysis and transient response analysis but must be 


included. A value of O. should be entered.” 
The only true variation in data record sequence between 
Small vibration analysis and transient response analysis 
comes in the next group of cards. The proper sequencing for 
ems Viorations 1S D-2, D-3, F-1; while for transient anal- 
ysis it is D-l, E-1, E-2, F-1. These are the Strategy 
Parameters and the difference lies in the need for integra- 
tion over time in the transient model. The D-1l record sets 
tolerances, limits the amount of CPU time, allows restarts, 


ween, while the & records control the actual time integra- 


; . 1s 
fen itself. 


Sein S Muse Oe a real number. Slee nts Ofvlen DplLind. to 


aeeeeeOt Submitted in the proper format. It will not record 
peer ror, but will Simply pass over the improperly input 
meee until finding Some in the format it seeks. This can 
Meer cO Sometimes dangerous distortions of the input 
records. 

Se tie wie Pate ties StUMdYy was Initiated it was 
Merete astvOood that implicit integration using K. C. Park's 
Permula (IMPL = 0, METHOD = 4 on the E-2 record) was the 
Only form of integration that should be used with USA. As 
of this writing, however, Lockheed informs that USA may now 
also be used with explicit integration. It has been 
Suggested that this might provide a faster and more econom- 
meal Solution. Due to time limitations that has not been 
tested here. 
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The D-2 and D-3 records are specific to modal analysis 
and deal with the number and range of eigenvalues to be 
Calendar oar 

Data records F through P pertain to the model itself. 
The records F through L provide information seneral toneae 
entire model and are listed in normal alphabetical sequence 
once through only, although specific cards may be repeated 
within the sequence. M through P, plus Q and R are specific 
to the individual shell and beam components and are repeated 
in sequence once for each component (See examples in the 
appendix). 

The F-1 and F-2 records provide grid point and stiffener 
Summaries, respectively. The information in F-1 controls 
generation of the finite element mesh. Inputs are Simm 
the number of nodal rows and columns for each shell.  AGiaam 
choice of which shell direction shall Gonsist of rows ame 
which columns iS unimportant as long as that choice remains 
consistent for the local X and Y coordinates across the 
entire model (rows are lines of consistent X, columns are 
lines of consistent Y). Entry is made sequentially; number 
of rows(I), number of columns( 1) for each shell ae orderaes 
input. The mesh created forms the master elements which 
will then be fitted into the structural cseometry by the 
information contained in the M records. 

The F-2 record is a summary of shell stiffeners. Use of 
this record could be extremely useful in determining 
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generalized responses on large scale models, but was incon- 
Sistent with the purposes of this study because of the 
inherent assumption that stiffeners do not deform. 

PeOorte esta DlLi Shes Cheweompatibdility conditions for 
mm@emevarious Shells. The instructions listed in the 
PeeiiC—-STAGS user's manual are straightforward, but refer- 
ence to Apvendix B may prove helpful for establishing proper 
linkages. It will be noted that in each case, for a flat 
plate, side 1 is that nearest and parallel to the global Y 
axis while side 4 has a similar correspondence to the global 
Peeeaxis . 

moceards Ia), hk, and b constitute the Data Table section 
of STAGS. Again the user's manual is self-explanatory and 
meell not be explicated except to note that the Material 
Record, I-1; the Material Property Record, I-2; and a 
Plasticity Record I-3 (if used); are listed in sequence 
Once for each of the materials enumerated in the B-3 record. 

If plastic behavior is expected in the model and the 
White-Besseling theory of plasticity is chosen (IPLST 
options 0 or 1 on the I-1 record) especial note should be 
made of the user's manual's warnings concerning the slope of 
meme ame COnnecring poants submitted for the discretized 
Stress-strain curve in the [-3 record. The change in slope 
at each point must be negative. Another point not mentioned 
in the user's manual, but which should definitely be noted, 
Meomeenatetie LTirsSt,» stress-strain data entry on the I-3 record 
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must be consistent with the elastic modulus entered on tne 
T-2 record. That is, the slope Of A2Whive Gravy wee aK 
Origin to that first data entry Should be Theva 
modulus. Also, while SpeGity ine ola stwe wie tes toe 
may actually decrease the run time for the TIMINT module by 
reducing the number of iterations required for each time 
step, it can increase the required computer space by more 
than 100%.7/ 
Tf plasticity is to be included in the model, there eae 
also some changes that should be made in the K records. 
STAGS determines the plastic responses of the shell by 
integrating through the thickness. The input NOT eee 
K-2 record establishes the number of integration points 
through each material layer. Since Simpson's rule is used, 
this number must be odd, NLIP may [te weer een yo ne mee 
(inclusive), but some care must be exercised in the choice, 


Since the larger numbers will add to the length of the 


an normal run for a model with about 300 nodes seems to 


be on the order of about two to two-and-a-half hours of CPU 
time for strictly elastic behavior. If the elastic range 
exceeded but plastic behavior is not specified, these runs 
can easily continue for 15 to 20 hours as the code attenome 
repeated iterations for each time step (the user can control 
this somewhat by specifying the maximum allowable time for 
any iteration), and the full number of specified time steps 
may not be attainable. Inclusion of Dlastaelt eee ee 
reduce the run time to roughly three to four hours. Opa 
of stresses beyond the yie'd point of the material, howeyveoume 
requires much larger solution data 1On022 2) Se  os 
both the TIMINT and POSTP runs. Typically, these have been 
on the order of 45,000 and 35;000 blocks Gestee ume slo 
the pressure levels used here. 
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MeN’ run; which, as has already been mentioned, is not 
Meebeniticant fOr the plastic case. The default value for 
Mets teley bse oe (his may be acceptable for models where 
the membrane stress governs, but will not be so if there is 
marked bending stress. 

Another change that should be made for runs including 
Siete tic responses iS to enter a value of i or 2 for the LSO 
entry of the K-2 record. Piste SswLSwiies COne thessA-SslAGs 
code will successfully complete an entire run, but the post- 
meee ssor, POSTP, will not print stress or strain ten 

The sequence of records M through R defines the proper- 
ties and loads on the shell units. As is noted most explic- 
itly in the user's manual, these records are read 
sequentially for each unit. All data in the sequence will 
be read before the M-1 record for the next unit is read. 

The individual shell geometries and the defining parame- 
ters by which they are to be connected to other shells (that 
is the means by which they are to be oriented to the global 
geometry) are specified in record M-l. EGEOB Ee opt lon. 0. 1s 
Meetul if the local coordinates of the unit already corre- 


meee precisely to the global coordinates. Ontion 1 is 


eee this may not be entirely true. iba OT dle raii.© 
get any useful response whatsoever out of POSTP, it has been 
Meeeccary to “trick™ the code into producing the desired 
stresses and strains by running the module with an 
menmooOZ. DAT fihke created by a STAGSI run for the elastic 
case. This has been true whether or not LSO has been 
flagged, with no noticable difference in the results. 


oe 


Similarly intended to be a labor saving device. <A caveat 
should be observed, however, that option 1 should only be 
used if the side connecting two shells is uniquely defined 
(i.e., no straight lines). It is strongly recommended that 
options 2 or 3, which allow more user control by Speci 7a 
boundaries through specific global reference points, 56 ee 
ine preilerenceswer o> eleven 

An option not iisted in the user's manual, Duly von 
could conceivably be useful for construction of Cylinder ae 
models is option 4. This allows the user to specify the 
location of a unit origin in the global coordinate 2 ee 
and then define the boundary in terms of one translations 
One rotation about the reference point. 

It should also be noted that STAGS will permit intveq@ias 
tion of shell geometries not included in the eleven which 
are indigenous to the program. This is done through a user 
written subroutine which may be specified by choosing ISHELL 
option 1 on the M-1 record. 

If the ISHELL = 1 option is chosen, the sheil geomet ae 
proverties normally input in record M-2 will be input) fie 
user-written subroutine LAME as stated above. Otherwise, 
Shell geometries are chosen from those provided by the code. 
Property values for these geometries should be input as 
specified in the figures of Appendax B. If plastic deforma 
tions are expected, particular attention must be paid to 
inputs to the M-5S record. First, simce Gatisstanmeuacdna cums. 
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is used to integrate across the mesh, STAGS will give erro- 
Meeus results (Or none at all) if that mesh is not located 
on the midplane of the shell thickness. ECZ must corre- 
eeemdingly be set to 0, and [PLAST set equal to 1. No 
further oes need be made on this record unless one of the 
Seer! units is located on a plane of Symmetry. eee nenpee S 
the case, STAGS will not compute lateral deformations 
Wienhout further information about how the shell is to 
feepond. ike ee ek ose Miauite ease onmuciod lia theses Oteminute 
randomly generated displacements or "random imperfections" 
eumene Shell nodes from the midplane. To trigger input of 
these imperfections, [RAMP should be set equal to 1, and the 
appropriate entries made on the M-6 and M-7A records. 

ie wser's manual's instructions and flow patterns for 
meerees N through R are clear and it is believed that, with 
Meeetole references to the input files in Appendix C, the user 
mm@emla have no difficulty completing the STAGS input deck. 

Me stated im the section on organization of the code, 
Peres Givided into four modules. Hach of these may be run 
independently once the needed data files have been produced 
Deememe previous module. Except for the required access to 
those generated data files, each module is run entirely in 
meama-alone fashion. Since modular input information is 
brief and error statements are generally good, this permits 
debugging without affecting earlier inputs and generates a 
eee eoeee OL Contidence in the modular output. 


V1 


It has already been mentioned that setting the IFLU 
Switch equal to 1 on the B-1 record of the STAGS1 run will 
senerate two data files FOROO3.DAT and FOROO4.DAT which are 
specific to USA. These are produced in addition to the 
FOROO2.DAT and output files which are normally generated by 
aol AGsiy runs 

Bach of the three data files FOROOZ.DAT, FOROCG y= 
and FOROO4.DAT must be renamed after the STAGS1 run for use 
in USA. The first, FOROO2.DAT, should simply be renamed for 
the file subdirectory containing the USA files (e.g., 
(*,USA)FOROOZ.DAT). FOROO3.DAT and FoORCC4 ius ie. 
renamed to correspond to the GRDNAM and MASNAM files speci- 
fied in the FLUMAS and AUGMAT modules, respectively. This 
may be done by entering a Simple RENAME command on the VAX. 
The user will then be queried regarding the name from whee 
and to which they are being changed. 

The GRDNAM file (the renamed FOROO3.DAT file) contains 
the global coordinates of the structural grid points genewe 
ated for the model by STAGSL1. The FLUMAS module then uses 
this information to construct a fluid mass matrix 7 ouiaee 
as well as generate the fluid mesh data (GEONAM) and trans- 
formation coordinates that relate the structural! ene 7 laa 
degrees of freedom on the wetted surface. The FLUNAM and 
GEONAM data files will then be used in the AUGMAT module to 


create a final augmented mass matrix. 
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ie Techs Cle) is phewcomponent of USA which assem- 
bles the data from the FLUMAS module and from STAGS1 into 
Mmeenwsoecitic constants and arrays required for the staggered 
solution process. The input deck for this module is very 
Short and experience indicates that ee With the module 
will generally be caused by input errors to either STAGSI1 or 
FLUMAS. 

Dea S time integrator module, TIMINT, links the results 
of the AUGMAT module which are stored in the data file 
Beye back tO STAGSZ and then unites with STAGS2Z to conduct 
a step-by-step numerical time integration of the governing 
equations for submerged structures exposed to shock waves. 
Wees time integration information input to TIMINT will over- 
See Carlier entries to the STAGS deck. It is suggested, 
however, that these entries should be made consistent. 

Input to TIMINT includes not only time integration 
Mmmioritacion but also the Spatial coordinates for the charge 
meet or the point on the structure closest to the charge 
(the standoff). These distances are made in reference to 
me@emoricgin Of the global structural coordinate as defined in 
Maemo lAGS input. One point of Caution: The surface-cutoftf 
provision in FLUMAS assumes a horizontal separation between 
MEMitemanad tarecer. As the Charge approaches a point verti- 
Cally above the target, some of the computed values will go 
Poesia pec avodie syne module to faal. The surtace-cutotrt 
mMinNCuLOn, GOeCS NOt automatically PiGieee Sei mOott and the 
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FLUMAS input must be changed accordamel, §ti a Vor wees 
aspéct is required. 

Provision has also been made in TIMINT to print nodal 
displacements without recourse Tema pOsmopreoee.ccr, It is 
recommended that this option not be exercised. A run for 
any of the other USA-STAGS modules can be made in a matter 
of minutes. Running TIMINT can easily consume several hours 
of CPU time for a moderately sized model. It is therefore 
not time or cost effective to run TIMINT each time a@ tev eee 
of displacements is needed. 

The solution data file, POSNAM, output by TIMINT is 
specific to the USA post-processor, POSTPR. If nodal 
displacements only are desired this is by far the easiest 
and most economical route vomrol lo] It is™not, howe v= 
the only option. USA-STAGS provides a selection of three 
possible post-processors. As mentioned, POSTPR 1s internal 
to USA and may be used to provide either data listings or 
graphical displays of nodal displacements. The other twee 
POSTP and STAPL, are a part of the STAGS package. A seG@eme 


SOlUtTIOnN Gata file Wines Lacs Pome” FOROZ2 .DAT, 18S crear 


aa is actually created by STAGS2 dureaeerie itis 


Ee 
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Pmeewmay be renamed to the STAGS subdirectory for use in 
either of THeselee 

Both POSTP and STAPL have been designed to provide not 
Only displacements but also stresses and strains; POSTP in 
mootlar tOorm and STAPL in graphical on a CALCOMP plotter. 
meee Lime Of This writing, the Naval Postgraduate School 
did not have a graphical capability on VAX/VMS so neither 
the POSTPR-graphical nor STAPL options were run at the 
Naval Postgraduate eenoolee = 

In accordance with common usage on the VAX, the sequence 
of commands required to execute any module have been entered 
into command files (*.COM). Examples are given in Appendix 
A. 

Once the appropriate file names have been entered into 
meemecOmmand file for the module to be run, the command file 
Memos Cxecuted by entering "@ command filename" for runs 
meammerne Lerminal; or “SUBMIT filename" for batch runs. A 
file definition statement (SET DEF [*.*]) is required for 
Memeot runs bul iS not necessary for runs from the terminal. 


205+ Should be noted that the TIMINT entry, NSAVER, is 


the frequency with which responses are saved for the USA 


meesporOocessor only. Ihe frequency with which data is saved 
ime ene POSTP and STAPL postprocessors is actually tied to 
mee restart data in NRESET. Pens tGesoe amar strain iwntorma= 


imeem tS required at each timestep, NRESET should therefore 
be set to l. 

-lpnosTpR Deol eons “SHOWN later Sseciions of this 
work were done by Lockheed. 


E.  DISCRETIZING thie ore 

As its name implies, the primary building unit in the 
STAGS code is the shell element. These elements may be 
constructed into shapes which are warped through almost any 
conceivable configuration in three-space and beam elements 
may be added, but the basic unit remains two-dimensional. 

The flat plate models constructed for this study, then, 
were ideally suited for discretization using Such Clem aes 
Each plate with its associated stiffeners was capable of 
being divided into several flat shell units composed of 
square or nearly square elements. Additionally, the even 
Spacing of the stiffeners allowed advantage to be taken of 
the plate symmetries to reduce the overall size of the model 
and corresponding length Of “tne ‘computes rum 

The transverse tee-stiffened plate established symmetry 
in two directions permitting a quarter model To seule 
That quarter model was then divided into five shells; two 
for the base plate, and three for the stiffener. These 
Shells were numbered beginning with the two shells 
comprising the base upward through the stiffener (Fig. 
4.2). This was done to permit changes to be made to the 
code input by simply adding or deleting shells. As long as 


the basic plate dimensions and symmetries remained the Same, 
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extensive rewriting could °C 4%Ot¢e¢ieu2el aa 
USA-STAGS modules.“ 

The STAGS code will generate the finite element mesh. 
The user need only specify the number of nodal rows and 
columns. In the case of the tee stiffened plate, an element 
0.6 inches by 0.6 inches was chosen in Order to keeo one 
mesh as fine as possible and yet be able to place the active 
filaments of the strain gages used during the underwater 
explosion testing entirely within the Pee eae This cores 
sponded to ll nodal rows and 6 nodal celinns sano 
unit 1, Li rows and 11 columns in Sheltie eo 
rows and 3 columns in each of the three stiffener shells 
(Pau re. 4.050 

Since the stiffener used on the second test plate was 
rectangular, three shell units sufficed to model it (i ace 
4.4). At the time that this plate stiffener was discre- 
tized, the limit set upon the number fluid mesh elements by 


the FLUMAS module was 171.74 The increased surface area 


oS Ome eee it was not possible To earry throu 


with this plan since the plate with the longitudinal reCtaue 
gular stiffener completely changed material dimensions and 
permitted only one axis of symmetry. 


23 


STAGS computes stresses and strains at the geometric 
centroids of the elements. Strain gages should therefore be 
centered on the element. 

24 


It turns out that this 1s relatively easy to change, 
but it involves technical considerations that are beyond ime 
scope of this work. QUeESUlOnS  2s0Ub Suen Motes eee 
Should be directed to Lockheed. The limit on most accounts 
has recently been increased to 400. 
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meecessitated by the use of a half-plate model Lasieeae. Of 2 
quarter-plate sev the total number of fluid elements beyond 
that number when the 0.6 inch element was used. Rectangular 
elements 1.2000 inches by 1.2857 inches were chosen for 

the base shells in order to allow strain gage placement at 
the center of each shell. Flements for the stiffener shell 
were 1.2857 inches by 0.6250 inch. This corresponded To 6 


nodal rows and 8 nodal columns for shell units 1 and 2, 


and 3 rows and 8 columns for oy cul Ween 6 anQldee ck ois: 





Figure 4.5 Half Model of Rectangular-Stiffened Plate 
For the final model of both the tee-stiffened and 

rectangular-stiffened patecw etme Pinive element mesh was 

paced at tne mid-planes of the Seles ee boo WAS OL 


el 


originally done. In an attempt to make the USA-STAGS 

Output correspond directly to the Strain Gage Outputs frem 
underwater testing, the plates were modelled with the meshes 
placed on the surfaces of the plates which corresponded to 
the surface of strain gage mounting. The half-thickness of 
the plate was then input as the ECG value on the M-5 record, 
thereby defining the true mid-plane for the code. This 
procedure worked extremely well for elastic runs. As stated 
in the previous section, however, the Gaussian quadrature 
that STAGS uses to calculate the plastic responses acroaag. 


the mesh surface is incapable of handling the integrations 


at any but the TEE 


For each of the test plates, maximum advantage has been 
taken of geometric symmetries in order to reduce the size of 


the finite element model. For the tee-stiffened plate this 


“oon ie important point iS reiterated here because it may 


not be immediately apparent from the STAGS run. The exper- 
ence gained here seems to indicate that when a negative 
eccentricity is input to the code, perfectly reasonable (but 
wrong!) strain approximations are obtained. The first 2aegeee 
cation of error in this met hogmrcn core ata ee 

and strain gage data arose when a run with a positive 
eccentricity was attempted. Under those conditions, STAGS 
produced no strain data at alls This Mineola, oie. eee 
not incapacitating, but it 1S inconvenient for 1t means gaee 
the data generated in the computer run will have to be 
converted to correspond to the strain gage information. ia 
theory this is simply a matter of multiplying the plate 
nhalf-thickness by the curvature Galewlated inst te seaweaae 
then adding or subtracting that trom the mid-plane Sera cae 
Given the massive size of the STAGS output, however, this 
means that either special processors will have to be 
employed or a conversion program will have to be written: 
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presents no further complications Since the resulting 
Quarter model is non-symmetric. For the plate witn the 
longitudinally oriented rectangular stiffener, however, the 
Mmesmit Of halt plate modelling 1S that the stiffener itself 
lies on a plane of symmetry. This presents a problem for 
STAGS. Although deformations along the symmetric plane will 
be calculated, no lateral displacements will be recognized 
Since, as far as STAGS is concerned, forces on either side 
are equal. This both affects the validity of the displace- 
ments that are calculated and causes a failure to indicate 
Stiffener tripping. STAGS solves this problem by allowing 
jmmemuser tO either input initial geometric imperfections in 
the shell reference surface or specify that the code 
generate random imperfections. AS stated in the previous 
section, this is done by flagging the IWIMP and IRAMP vari- 
ables on the M-5 record and then following the flow direc- 
fens Listed in the manual through the M-/A record. For the 
rectangularly stiffened plate used here, imperfections in 
the first three harmonics of 0.001 inch were eee 
morally wsalthouch met Strictly related to plate discre- 
ieee@acion, there iS one other point not mentioned in the 
STAGS user's manual, but which 1s essential if the computer 
mete 15 CO successfully compute plastic deformations. oe 


oces the STAGS1 inputs for the rectangular stiffened 


plate in Appendix C. 
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plasticity is to be included, the variable IPLS? on thew 
record must be set to some value between O and 4 to indicate 
which of four plastic Strain heen ees es geet ee iseceem a 
non-zero positive value must also be specified for the vari- 
able NESP. This NESP value specifies the number Of Gauge 
points on a stress-strain curve for the plate materia lain 
will be input to the I-3 record. There aré two points ae 
the user's manual fails to make clear about these stress- 
Strain points. First, this curve requires Stress ees 
neering strains beyond the yield point. Second, the first 
point of the curve must correspond closely to the moci Bites 
elasticity of the material specified in the [-2 record = 
All the material and stress-strain characteristics for age 
6061-T6 aluminum used to manufacture the test plates here 
were taken from the Military Standardization Handbook, 
Metallic Materials and Elements for Aerospace Vehicle 
Structures [Ref. 14: p. 3-184]. 


ane default value, O, automatically triggers use of 


the White-Besseling theory, as does the set value l. 

Tere 1s, the slope from the origin to the first se@auue 
specified must be the modulus of elasticity. The White- 
Besseling theory will apply approximations to attempt a 
patch-up. If the discrepancy is large, however, continuity 
between the elastic and plastie responses will be lost ae 
the resulting output will be in error. 
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Peewee MODAL TESTING 

Pepermination Of a Structure's mode shapes and natural 
frequencies can provide a first test of a finite element 
Meee Ss Ability to predict the dynamic responses of an 
excited system. Perhaps of even more importance to the 
purposes of this study, it can be an excellent check of the 
miewervural data input to the code. Failure to provide 
Proper connectivity between shell units or inaccuracies in 
mem anOuUtL AGimensions or material properties will show as 
sross discrepancies in the modal analysis. 

The Naval Postgraduate School is fortunate to be one of 
the few institutions to have a Hewlett-Packard System os451C 
Member Analyzer with modal testing capability. 

The technique used by the H-P 53451C iS in many ways 
Similar to the solution processes used by the finite element 
codes themselves for it involves solving basically the 


Same governing equation of motion [Ref. 15]. 
MY + Cy + Ky + f = 0 (eqn 5.1) 


Peetiee WOtLh ¥ and ~f are functions of time, the transfer 
mmmection, h(y), may be found by taking the Laplace trans- 


minor KQOUALION O.1 


eo 


B(s)X(€s) = FCs) (eqn 5.28 


where, 
2 r 
BCs) = ‘Neu eo (eqn 5 ae 
and, 
eo z 
H(s) = BCs) (eqn 5.4) 


If the exciting forees and displacements are Hie wie 
may be readily found by solving the equation which resuiilee 


from substituting Equation 9.4 into Equation oOf2 > temo 


H(Cs)F(s) = X(s) (eqn 3S.om 


For an oscillating system of order n, the transiter jie 
tion will have 2n poles occurring in complex conjugate pairs. 
Fach pair of poles will cause a mode of vibration fie. 
structure. Each of these poles iS a complex number whienl 


may be expressed as 


Pi = =o Fey ic (eqn 5.6) 


where - 1s the damping coefficient and Wy the natural 


frequency for that mode. 


Once oy and Wy have been determined, the resonant 


frequency 8, (in radians per second) and the damping ratio 


k 
Cc follow directly. 
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: > ak ae (eqn 3.7) 
v4 = JJ 7 ne a Aref 
i q 

oe 

me 


The exciting forces and displacements that must be 
provided to the system may be found in a number of ways. 
The method chosen here was to excite the olate using a 
modally tuned impulse heummese 7 © An een was 
attached to the plate to record the resulting accelerations, 
from which displacements were found by numerical 
ieee cration. 

The two fundamental assumptions of this test procedure 


(Ref. 15] are that: 


ieeevodgal frequency and damping are constants for any 
meansrter function taken from the structure. 


2. Modes of vibration can be excited from anywhere on 
an elastic structure except at their node lines (Where 
nMOme<ere ac 1Onewls DOSSIDLe) ; 

It should therefore be possible to describe all modes of 
Meemerest by Simply choosing a sufficient number of excita- 
tion points on the structure and by placing these points 
close enough together to minimize the possibility of 
Zee cation. 

For the two plates used in these tests (the rectangular- 
and tee-stiffened PEeStemMOdelS),= mine points of excitation 


Senecio ceos. 


i aceemodel S0oA07 . 


S7 


were chosen on the back of each eee Care was taken to 
avoid symmetric placements of the enters The accelerometer 
was fixed at a single location arbitrarily chosen but Spaeee 
away from the clamped boundaries. 

When any location was tapped with the modal hammer, the 
resulting forces were fed back through the hammer to the 
Fourier Analyzer; as were the accelerations from the accel- 
erometer. The Fourier Analyzer then computed the displace- 
ments and solved the Laplace transform to create the plot 
Shown in Figure 5.1. Each of the peaks in this representa- 
tion corresponds to the poles of the transfer funct 10m ig 
Single point of excitation. By then combining the DlOtSiiaas 
all voints, an accurate representation of the system modes 
was obtained. Tt was then simply a matter of analyZzino eee 
of the poles in descending order of magnitude to whatever 
total number of modes is of interest to the user to deter- 


mine the dominant modes of the system. 


B. UNDERWATER SHOCK TESTING 

Underwater explosion testing for all studies conducted 
in this series to date has been done at the West Coast Shock 
Test Facility (CWCSF), located at the old Bunter ae cea 
Naval Shipyard near San Francisco, Calivepiia) 9! aise 


Se olates were placed face downward upon a solid 


marble slab. 
32 
The actual number and Spacing Of The excita ton pom. 
seems to be primarily @ matGer of §o4e ewe 
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Navy activity reporting to the Supervisor for Shi pone awe 
Conversion, and Repair, San Francisco. As the only licen 
facility of its type on the West Coast, they have gained a 
considerable amount ot experience in conducting underwater 
explosion testing; although most of it has been large scale 
work in support of the Navy's shock qualiticgalion Sues ae 
rather than experimentation. 
1. The* Testo Place rewm 

In OFder tO PROV tae ar me mere backed condition thee 
was required to simulate a ship's hull and to ensure rigidly 
clamped boundary conditions, the test plates used ince 
series of studies were securely bolted to a heavy steel 
backing structure (Fig. 5.2). This structure Cconsisteqmeuaee 
box constructed of 14 inch thick AG structural steel ie 
wide flange was welded around the open side of the box to 
which the edge of the plate was bolted. An O-ring.-gasket 
was fit into a channel machined into the Suritace Of Sene 
flange to provide watertight integrity. Access £0 Insta 
mentation on the inside of the box was allowed by a subma- 
rine style penetrator affixed through the bottom.” 

To perform the actual testing, the plate and backing 
Structure were suspended by heavy steel cables below two 
large pneumatic floats in the configuration developed in 


ede a more complete description of this backins straee 


ture and the rationale used in designing it, see Reference 
3, pp. 88-94. 
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Chapter 2 (Fig. 3.3). The cables wereweut (6 feta 
assure that the nearest edges of tne test panel Stifteneme 
were maintained at depths as close to thirteen feet below 
the surface as possible. Once the entire apparatus had been 
towed into the bay, the explosive charge was lowered into 
position nine feet above the stiffeners in line YVeuur ween 
With the Center "Obese nemo la tee 
2.° LS tains: 2 red 
In his work, which has already been cited several 
times, LT Rentz outlined the basic theory governing the test 
of the flat plate and began developing the procedures which 
will be enlarged upon here. 
Essentially, the information sought for coOmparuecs 
with the computer simulation codes being tested has been: 
1. Fluid pressures at the plate and in the freefield. 
25 oOotvrain 17 @eme Plates 
3. Final deflections ln Sino we a2 
Instrumentation used to gather this data was as follows: 
a. Pressure Measurements 
Three vressure gages were used in the tests 


conducted here. Two of these were located just as they had 





been in the previous study [{Ref. 3: p. 46]. One was located 
one foot above the center of the plate and the other was | 
clamped to the aluminum block bolted to the side of the 


plate. The forward edge of this second pressure gage was 
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positioned so as to be on an are of the same radius as the 
standoff distance. 

The pressure recordings that had been obtained 
from these two gages had been so subject to interference 
from the vlate structure and the rigging of the ves so ae 
form that the futility of conducting further tests) wittogeee 
true free field pressure gage was readily obvious. Accor- 
dingly, although the first two gages were retained as 
backups, a third gage was positioned on an arc equal tothe 
Standoff distance but removed five feet from the centerline 
of the plate (Fig. 5.4). As will be discussed in the 
following chapter, the resultS obtained from thes y2ase 
during the test shot on the tee-stiffened plate were 
excellent. 

All #@sts in this Series mtr) meee one-quarter 
inch tourmaline crystal pressure gages rated to 10,0003. 


The response Eon RO? 1s 42 fimeticn GCiuene Geen: 


constant for the charge, 3, divided by the transit @ imewen 
the shock wave across the gage, ty [Ref. 16]. For the 
eight pound TNT charge chosen as the standard for these 


studies, 


9 = 0.15 msec 


= ermine as the ratio of the apparent peak pressure to 
the actual peak pressure. 
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From the chart on page 17 of Reference 16. 
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This indicates that a very good correlation can be expected 
between the actual and measured pressures, and is also the 
best indication that the correct gage size is beins Wee 
b. Strain Measurements 
Perhaps the single most critical aspect of 
Obtaining successful data from underwater explosion testing 
of the type conducted here is the choice and application of 


the strain gages. 
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For the tee-stiifened plate, twenty SR-4 oe 
FAER-25-12-S13ET strain gages 2 were attached as snown in 
Figure 5.5. These were the same type used in Reference 3; 
memes | tCn@uech rated to only 2000 microstrain, it had 
been demonstrated that their range could be extended suffi- 
ciently to provide useable data over the very earliest 
increments of the shock response. As in the earlier test, 
BLH's SR-4 Permabond 910 was used as the fixative, and the 
attached Strain gages and their wiring was covered with 
PR-1422 A-1 Base Compound, manufactured by Products Research 
mmeectiemical Corporation of Glouster City, ee This base 
compound was also used on the plate with the rectangular 
stiffener where 19 MX type PAHE-250BA-350EN strain eee 
meme. 3.6) rated to approximately 200,000 microstrain were 
Used in an attempt to increase the time span over which 
information could be gathered. Also, at the recommendation 
epee de GAistributor, BLH type EPY-150 epoxy adhesive was 


used instead of the Bondo 910.°° 


Soienuactured Syeetiietbecet DOonmLCS of Waltham, MA. 

Tae was done to provide and smooth compliant 
covering that would keep the gages and wiring from being 
torn off with the passage of the shock wave, and to seal the 
gages from saltwater once the rig was submerged. 


"Manufactured oy eoeeemimencal otvress Analysis Products. 


on theory, the epoxy provides a stronger bond than 


does the Bondo, which 1S a contact-type adhesive. 
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In practice, although some useful data was 
obtained from each test, neither was a particule 
Successful source of Strain initormat@on. The numoewae. 
gages lost far exceeded the number thal, perir6c71 a. 

In recapitulation, a number of fault ee 
found for which the following corrective measures have been 


proposed. 





1. Although the MX type strain gages greatly exceed 
the strain requirements expected in these tests, 
their use is strongly recommended. Little or no 
loss of accuracy 1S experienced in the lower strain 
regions through use of the high-elongation gages, 
whereas mid-range gages tend to peak right at the 
Strain ranges expected. 

2. On any Strain gage, the weakest points are the sold- 
ering tabs to which the gage wires are attached. 
These are brittle and will break if any bending 
stress 1s applied. All motion should be restricted 
to the plane of the gage. This means that the strain 
gage wires must be extremely flexible so that any 
bending motion, either during installavlon ere wee 


result of plate deflections, is not transmitted to 


the gage tabs.  Ratedwstrain Gace (aera 
recommended. ~~ 
on the tee-stiffened plate, 28 gage wire made specifi- 
cally for strain gages was used. This prevodi eee ee autre 


LOO 


3. Even with flexible strain gage wire, a determined 
force along the length of the wire, as would ve experi- 
enced during plate deformation will break the tab, 
Pheu wiCe wOrmm Ene sSOlucderme@emnectlon between the two. 
Although a plastic medium, the PRC Base Compound 


used to protect the gages and wires was not flexible 


enough to allow the wires to expand with the plate. 7° 


A suggested solution to this problem is to lay the 
gage wires with some s-curves to allow elongation as 
the plate deforms. The wires should not, of course, 
be otherwise affixed. A piping of Barrier WD, a 


tacky, semi-fluid substance manufactured by BLH which 


Satisfactory as far as the gages themselves were concerned, 
but presented definite problems on the human side of the 
eGuacion. The wire was so fine that it was impossible to 
orem With it. It also proved to be highly vulnerable to 
damage while the test structure was being transported or 
mememeod. For the following shot, the attempt was made to 
use 20 gage bell wire. This withstood handling magnifi- 
cently, but was too stiff for the soldering tabs. Annealed 
42 Or 24 gage wire designed specifically for use with strain 
Mees appears to be the best possible choice. 
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This base compound Seems to be an extremely protective 
feeeeim. A couple of cautions should be issued governing its 
use, however. First, while the compound must be used within 


One hour of mixing base and accelerator, enough of that time 
Should be dedicated to stirring to ensure that the compound 
memeenoroughliy mixed. Although this observation is not 
Eomoeletecly confirmed, it appears that insufficient mixing 
leaves a residual which can permeate wiring insulation and 


actually destroy very fine copper wires. Second, the 
compound 1S not compatible with the Gage Cote lacquer 
normally painted over strain gages as a sealant. If the 


Mes -wcompeund 1S TO be applied directly to the strain gages, 
Gage Cote should not be used. 


LO 


looks somewhat like a combination of beeswax and 
petroleum jelly, could then be eis over both gages 
and oe This piping of Barrier WD should then 
covered py a layer of aluminum fol1l and then the whole 
overlayed by the usual PBC Base Compound. The net 
effect is an armored "tube" of Barrier WD through which 
the wires can run when the plate deforms. 

4. The fixative used to attach the strain gages should, 
itself, be a subject for careful consideration. Even 
though the epoxy is supposed to be the stronger medium, 
experience gained here indicates that it 1s not as 
Satisfactory as the Bondo for explosion testing, a 
postulated reason for this seeming discrepancy 1s that 
the epoxy places a thicker layer of adhesive between 
gage and plate material. Given the extremely rapid 
loadings experienced, intense shears are created in the 
thicker substance and the epoxy literally cracks and 
spalls off. The Bondo 910 therefore seems to be the 


preferable Cee 


7 bakers pastry tube works very nicely [0 dour 


*2No failures due to loading were noted wii he). a aosme 


gages attached with Bondo. Therefore, Even if the cracking 
problem had not occurred in the epoxy, the greater separation 
between gage and plate material that the epoxy creates would 
argue in favoew ot the Bonde: 
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c. Plate Deflections 

Plate deflections were by far the most reliable 
Source of data that could be obtained as a result of testing 
Since no electronics were involved. To get these, precise 
measurements were made at specified points on the plate both 
before and after the test shot. A jig was devised which 
took advantage of the rigid plate boundaries to provide a 
Mmepeerence table for a dial indicator depth gage (Fig. 5.7). 
All measurements were made while the test plate was bolted 
Onto the backing structure to ensure that plate warpage did 


not affect the readings. *° 


See DATE RETRIEVAL 

Signals from the gages mounted on and about the test 
plate were fed back through shielded cable to a common 
floating terminal box. From there, they were passed to an 
iistrumentation shack located on one of the jetties defining 
the small bay in which the testing was conducted. A sche- 
Merc Of the setup iS shown in Figure 5.8. Calibration 
levels were applied to all pressure signals by Endevco 
Meee! 4470 Signal Conditioners. These same signal condi- 


tioners were used to provide calibration for all but twelve 


43 yeflections have already been given for the rectan- 
Mmumetr Stittened plate in Figures 3.7 and 3.8. No post-shot 
measurements were made of the tee-stiffened plate because of 
the heavy damage which was inflicted when the plate was hit 
Beeone Of the Supporting cables. | 
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Figure 5.7 Jig Used to Measure Plave Berlecr1 on. 
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of the strain gage channels. Those remaining twelve 
channels were calibrated by two six-channel Honevwell ‘iodewl 
B2-6 Universal Bridge Balance units. The strain gage signal 
were then amplified by Bell and Howell CEC 1-168 DC ampli- 
fiers before final recording. yee wi ees were recorded 
on two twelve channel Ampex FR 1300 tape recorders. 
Calibration levels for the tourmaline crystal presse 
sages were all set at the 10,000 psi theoretical maxima 
For the strain gages, however, attempts were made to use the 
computer codes to provide estimates Of expected sca 
uSing the pressure and decay Constant GCaleulateqd ia 
Equations 42.2 and 2.3. This Was €svecially crit ten eee 
the strain levels for the tee-stiffened plate where some of 


the strain gages were expected to be operating at or beyond 


tigear rated limats. 


Unfortunately, difficulties encountered in learaine 
to use the USA-STAGS code prevented runs from being-made 
before the shot dates.*" 

Advantage could have been taken of earlier estimates 
used for the plate tested in Reference 3 if the charge 
weights had remained the same; however, such was not the 
case. 

Since 8 pound TNT charges are listed as Navy standard 
stock items and order had been placed well in advance, no 


44 
In fact no successful runs were made until lone aftes 


all experimentation had been completed. 
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emeolemsS were anticipated in odtaining charges of this size. 
As the date of the first test neared, however, it became 
obvious that no 8 pound charges would be available. The 
nearest alternate that could be provided in time for the 
scheduled test was a 94 pound TNT charge. 

It can perhaps be appreciated that dimensional changes 
Mempne Lest Dlatfiorm are not simply a matter of moving the 
charge location. As charge weights are changed, the stand- 
off distance and the charge depth both must be changed; 
fimcheanvolves cutting new supporting cables and recali- 
brating flotation depths. All of which is extremely time 
consuming. When the personnel at WCSF were contacted about 
Ig@en possibility of changing the configuration for a 93 
memma Chargce, they advised that it could not be done in time 
for the scheduled tests. 

The choice for the first test was clear, either run the 
test with the existing platform and a 93 pound charge, or 
abort the experiment. It was determined that more benefit 
@eomee be gained from conducting the test. 

Wide values for Po and 9 found from Equations 2.2 and 


QO 


Emo tor a 94 pound charge are: 


Po = 4081.7 psi, 9 = 0.161 msec 


ieee Peep eessure was markedly higher than that for 


mies pPOUNd GCharcge. The maximum bubble radius, too, was 


LO7 


much larger. ren boQuamtene 2 ee = 14.03 ft. This 


gives a charge depth to bubble radius ratio of 


4 ft 


oft 285 
Teese 


which implies that some pressure release will take place 
into the atmestomere. 

For the second test shot, the problem was quite 
different. As was discussed in Chapter III, the tee- 
stiffened model appeared too rigid to successfully demon- 
Strate stiffener tripping. The plate with ane longit ud imma 
rectangular stiffener was a largely tentative attempt to 
produce that action. It was desired to view the responses 
of the stiffener under a variety of loadings, beginning with 
small charges and concluding with ones in the 8 to 9S pone 
range. Again use was made of what was available. in thao 
case, the available charge weight was 4 pound. 

For a = pound charge, the maximum bubble radius is 3.02 
feet; which leads to a charge depth to bubble radius ratio 
of 1.33. This implies that bubble pulsation wm yew 
problem. Fortunately, however, the time of the pulse is 
relatively late, t = 0.17 seconds. Thus, by Pimtetnse rae 
period of study to early times the bubble pulse can be 
effectively neglected. 

Calibration levels for the test shot on the tee- 
stiffened plate were determined by making an EPSA run on the 
plate model used in Reference 3 but with a 943 pound charge 


ers 





weight. As has already been stated, geometric similarities 
between these two plates permitted broad generalizations to 
be made between them even if exact comparisons could not be 
Seew. Che results of that run are shown in Table III. 

No similar comparisons could be made with the plate 
with the longitudinal rectangular stiffener. Strain gage 
calibrations were therefore arbitrarily set at a value at 
the higher end of the strain gage range which could be 
expected to exceed all requirements. In this case, that 
value was 150,000 microstrain. 


The results of these tests will be discussed in the 


next chapter. 


Table I1{I 


Calibration Data for Tee-Stiffened Plate 
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Gage BPSaroerain* Gage EPSA Strain* 
i 62 e Hea 76800 
2 2n® 2 7610 
> Toes Ike 7610 
4 8500 14 Soke 
2 ae? OO ne 7610 
6 69400 16 Zi Oe 
i 69400 17 6210 
8 69400 18 6210 
9 LOOOO0O ie) aoe 

FO Loo o 20 S5200 


aeeeevalues In mMicrostrain. 
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VI. RESULTS “AND CONGEe tert: 


A. MODAL COMPARISONS 

As stated in chapter four, the primary purpose of modal 
testing is to provide a check of the finite element model. 
A properly constructed model should exhibit substantia 
the same fundamental frequency as that determined by modal 
testing and should demonstrate close coincidence with oases 
frequency values across the range tested. Exact concurrence 
cannot be expected since the computer code 1s analyZino@eane 
idealized model which will not be affected by the externali- 
ties inherent in the actual test. Also, the numser ven 
frequencies determined by the computer is a function of the 
number of iterations than can be made over a given time 
Span. Certain economies in the use of computer time must 
therefore be observed. While the code might be able to 
calculate the frequencies for a large number of modes; such 
calculation will be purchased with a heavy expenditure in 
CPU time which 1s not really necessary for an adequate 
check. 

Accordingly, the scheme followed here has been to limit 
the number of modal frequencies computed by the code to 
those that could be determined anmewoemuns [en eacnee aaa 


Kach of these runs can be performed entirely withiaeine 
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combined STAGS1/STAGS27° processor Without reference to the 
Meeadd-on. The first run Simply specifies a broad range 
Within which the computer can search for eigenvalues; the 
second uses one of the mid-range frequencies determined from 
Meewetr St run Gigenvalues as a ‘seed” for further iterations 
in order to produce a refined product. 7° The actual perform- 
ance of these functions is controlled by the D-2 and D-3 
mecords ofthe STAGS input deck. Since the procedure is 
outlined clearly in Reference 13 and example input files are 
provided for both the tee-stiffened and rectangular stiff- 
ened plates in Appendix C, further elaboration will not be 
provided here. 

Modal testing of the milled plates was conducted in the 
mammer described in Section D of Chapter IV. The results ofr 
both that testing and the STAGS2 estimates are provided for 
comparison in Tables IV and V tee-stiffened and rectangular- 
Stiffened plates, respectively. 

It will be noted that, for each of these plates, the 


STAGS2 predicted values seem To correspond very well with 


“Generally referred to here as a STAGS2 run. 

eee Ms els nod trestinwwonr the milled 
plates was conducted prior to the computer runs, essentially 
the same frequency range was used for the computer runs as 
Mere pme testing. Over this range, the first STAGS2 run will 
compute five modes with their corresponding natural frequen- 
aa Lhneousn the uSer's Manual States that calculation of 
up to twenty eigenvalues can be specified, and that all 
Sxisting eigenvalues up to that number will be computed, 
five seems to be the true maximum. No increase in the 


ia 


Table IV 


Modal Comparisons for Tee-Stiffened Plate 


STAGS2 Predicted 
Modal Frequency (HZ) 


Measured Resonant 
Frequency (Hz) 
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Table V 


Modal Comparisons for Rectangular-Stiffened Plate 


STAGS2 Predicted 
Modal Frequecy (Hz) 


Measured Resonant 
Frequency (Hz) 
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values from the modal testing, tne greatest relative error 
(using the modal test value as a base) being 6.32% for the 
tee-stiffened plate and 2.16% for the rectangular-stiffened 
plate, and the average relative errors being 3-16, and ie 
for the tee- and rectangular-stiffened plates, respectively. 
Less promising is the seeming failure of the code to 
predict the fundamental frequency for either of the plates. 
It is suggested, however, that this is not due to any errors 
in-the code or in the inputS to the COde veers.) 1c meee 
attributable to the fact that the DounGdar, Beemer en 
actually tested were not the clamped-clamped conditions 
specified in the STAGS inputs. The H-P S451 yieu. 
Analyzer appears to be quite sensitive to relatively 
small changes in the boundary conditions. When the tests 
were conducted, the plates were laid face-downward upon a 
heavy marble slab. No attempt was made to further clamp 


the heavy aluminum frame into which the test plates are 


amount of CPU time allowed for eigenvalue analysis (NSEC on 
the D-2 record) seems to affect that number. Computatiommes 
elgenvalues over a range without any seeding of Known values 
1s a slow process prone to convergence problems. It showed 
therefore only be used to derive a value which can be used 
for the second run; or it can be eliminated altocetner ae 
one of the natural frequencies is already known. If sucha 
frequency has been determined, it can be entered directly 
into the code and all frequencies closest to muitiples cme 
fixed interval (SHIFT on the D-3 record) will be computed. 
This process is reasonably fast and will produce ail 
existing frequencies up to the number (NEIG) specified. 
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Goeed, nor was it possible to ensure an absolutely flush 
mount between plate and slab. 


When the tee-stiifened plate was tested, the bolt noles 


had already been drilled in the frame. These bolt holes 
seemed to cause reflections which were picked-up by the 
Fourier Analyzer as additional modes. These may be observed 
to some degree in Table IV in the repeated tendency for a 
group of frequencies to cluster about a single value. The 
effect was far more noticeable on the oscilloscope traces 
Sereinea aS GCach Of the excitation points on the plate was 
leced. When all of the variables are considered, it is 
believed that the actual correspondence between the modal 
test results and the STAGS2 prediction is quite good and 
Moraes adequate confirmation of both the code's ability 
memoredict the natural frequencies of a simple structure 
memeort the finite element models input to the code for this 


Seudy. 


Peo OMPARISONS WITH UNDERWATER TEST RESULTS 

As has already been stated in Chapter V, all underwater 
Pi@ee. teSting for this series of studies was conducted at 
Mmeemyvest Coast Shock Test Facility. The outputs from the 


Si is strongly recommended that future modal testing 


be conducted before the bolt holes are drilled in the frame. 
An attempt might also be made to weight the frame with sand- 
moo FO Simulate aS mearly as possible the clamped boundary 
meer tOns input to the code. 


ale 


pressure and strain gages used for the tests along with 
their associated calibration and time signals were recorded 
by two Ampex FR 1300 tape recorders. After each of the test 
shots, the data tapes and one of the tape recorders were 
transported to the Naval Postgraduate School where the tapes 
were played-back, filtered, and recorded on disk on the H-P 
9451C Fourier Analyzer for later analysis and pleceunee 

As thorough a pretest analysis as was possible without 
actual disassembly was conducted upon the equipment at the 
West Coast Shock Test Facility before each shot. Calibra- 
tion levelseiorra Senene devices were also checked 
using technicians and equipment provided by the Naval 
Postgraduate School. Despite these precautions, the Quai 
of data obtained from the tests was highly problematical. 

An axcerieet free-field pressure history was obtained 
from the test shot conducted on the tee-stiffened plate 
(Fig. 6.1); fully compensating the difficulties encounteme 
in mounting the third pressure gage in a position where it 
would not be affected by interference. AS can readily be 
observed, this curve displays all of the classical expenmeae 
tial decay characteristics predicted in Chapter Ll. Site 
this is so, the decay constant will correspond apoproximageue 
to the time it takes fOr (Che Dressite tesco o sr. = (abou 
.36788) of itsS peak value. Prom Preure 6. ee pee 
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Freefield Pressure for Tee-Stiffened Plate 


tose 


(3125 psil)(.369/7638) >= 136 eee 


which occurS at about L706 mii -eeae 

These values were input to the USA-STAGS code to conduct 
a post-shot analysis for comparison (060 (Test Gaca- 

The peak pressure and decay constant obtained from the 
test results are much smaller than those predicted by 
Equations 2.2 and 2.3. In retrospect, this 1S note sae 
cularly surprising. First, it will be recalled that Equa... 
2.2 and 2.3 are derived from empirical resuitsS and are Ga 
fore subject to wide variations. Even more pertinent, it 
was found in Chapter V that the ratio of the charge dep aiuaae 
the maximum bubble radius was much less than the .50 
threshold for pressure loss. The decreased pressure and 
more gradual decay are exactly what would be expected from a 
charge that had vented a Significant portion Of i1tS eee 
in the first moments after detonation. 

Given the good results of the first shot, a high ¢eawe 
of optimism was established for the results of the pressure 
data from the second, the $ pound charge on the plate Wilag 
the longitudinal rectangular stiffener. That optimism se 
not JUStifweds 

Upon playback of the data for the second shot, it became 
clear that the amplifiers from all three of the pressure 
Spaces had become totally saturated. This was surprising 


Since the pressure gages were calibrated for 10,000 psi and 


DS 





lh 


the peak pressure for a pOouUnG@e IND charge is far iess 
mean that. 
The methods by which the two tests were conducted 
Suggested a probable cause for the pressure gage saturation. 
Once all wiring was connected and the gages tested, the 
meu plattorm was lowered into the waters of the bay. Heavy 
Pemeyeards were attached to either end of the floats, one to 
be used aS an outhaul and the other as an inhaul. Using 
these, the platform was positioned in the bay a safe 
distance from the jetties before the explosive charge was 


installed and the test conducted. 


Me conduct the test with the 9 


bole 


pound charge, the 
platform was centered in the bay to keep it as far away from 
the jetties as possible. There was a heavy wind blowing and 
when the cable used to carry the flring current to tne 
explosive charge was taken re to the platform it was blown 
well away from the platform and the instrument cables. 
Since it was separated from the other instruments, the free- 
field pressure gage cable was farthest of ail from the 
detonating cable. 

When the test shot was conducted on the rectangular- 
Stiffened plate, the charge was sufficiently small that the 
platform was kept much nearer plierside to make it outa cx © 


menaie. All of the cabling was run out to the platform 


erene the inhaul. 
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The voltage used to detonate the explosive charge is 
quite large (800 volts DC) and sets up a noticeable induc- 
tance on all fe eewnin dee eee Ee VAC an ripe In tac 
on the strain gages, which are very low imoedance devices 
and therefore relatively insusceptible to inductive noise 
this was used aS a convenient, 11 not Wery AaeCuragoe- ae 
time nea 

The pressure gages, on the other hand, are extremely 
high capacitance, high impedance gages and correspondingly 
are very much affected SaaneaRe currents. It 1s thereiewe 
believed that the close proximity of the detonation cable 
was the sole cause of the saturation observed on the pres- 
Sure gages during the second test. Be that as it may, no 
pressure history was obtained for the rectangular-stiffened 
plate, making it necessary to use the peak pressure and 
decay constant calculated from Bouations 22 ade 2. oan 
the subsequent USA-STAGS run. 

The lack of a zero time Signal was not the only Gera 
ciency noted on the strain gage histories. Three others 
Should also be mentioned by way of warning for those 
CONdCUCT ING asim lar Sees 


48nhis was a major deficiency in the way testing was 


conducted at WCSF. No accurate zero time Was DrovrTeca om 
any channel. The detonation signal was used as a best 
alternative, but Since that typically @eover-cea eo ee 
about one millisecond it was almost useless when only the 
first few milliseconds were being studied. 


eZ 


Before the test platform was placed in the water, all of 


lee 6 ene le cumanenel 


ct 


the strain gage signals were balanced with 
resistors, setting their Signals to zero. AS soon as the 
platform was placed in the water, the change in temperature 
effected a corresponding change in the resistance of the 
strain Sages. (his in turn cnanged the balance of the 
bridge. This was an ongoing problem throughout the period 
of the tests; many of the Strain gages never were brought 
[eporproper balance. The problem was exacerbated by the 
fact that the strain gage and its associated wiring formed 
an extremely long leg of the balance bridge which was 
Subject to the environment of the bay. The balancing leg 
Semstisted of a calibration resistor located in the instru- 
Memeestiack; a very different environment from that outdoors. 
The second of these three problems 1S one which is 
readily ees by the experimenter if he is aware of 
it. When the test plate and backing structure were sent to 
the West Coast Shock Test Facility, all strain gage leads 
were connected to junctions in the terminal box located on 
m@emsicde Of the backing Structure. Connecting these leads 
through long cable runs was one of the many tasks that was 
performed by WCSF personnel between the time that the plate 
was delivered and the date of See cae Shot. Their experli- 
ence was primarily with the Navy's shock test program, 
wherein the magnitudes of the strain outputs were of more 
Value than the sign. Therefore, little attention was paid 


ab 


to the polarity of the strain gage connections. Thus, the 
histories recorded can be the exact inverse of the true his- 
tory. This can bé especially debilitating when comparised 
of the shapes of strain curves 1s of particular importamee 

Even more debilitating, but far less obvious, was the 
last of the problems noted. 

The amplifiers on the Ampex FR 1300 tape recorders used 
to record the strain signals are calibrated for a 40), doeugeee 
tion from a one volt rms input signal; peak-to-peak. a7 
maximum signal deflection that 1s theoretically possibi-w 
about 1.4 volts in either the positive or negative direction. 
The calibration levels input to the system are based upon 
this maximum; Signals greater than 1.4 volts will Satu@awe 
the amplifier and the playback will simply be a straight 
line at the maxima. Yet on the recordings for vinguaie 
all of the strain gages which Survived Geom eueogmee litt =. 
Wiese records, the output was greater than the 1.4 volt 
limit, but no concurrent indication of amplifier SauUgauameee 
Figures 6.2 and 6.3 are offered as typical representations 
of this phenomena for eh Oieb heat wemekatves: 

A number of possible sources for the strain gage error 
were considered: 

ls The Gai bration Sine were improperly set. 


2. The signals were being altered either as they were 
being input to or owtpuc from tne Fourter 2na) 2a 


3. The playback amplifiers on the tape recorders were 
OUT Of Ca blorawe ne 


dee 








Figure 6.2 Untranslated Strain Output for Tee-Stiffened 
Plate (Gage 11). 





Fueunpe 625) Unt Pama eainee Strain Output, Rectanzulams 
Stiffened Plate (Gage 15). 
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4. The calibration signal was being -input to the signal 
CoMUmerreners at une wrens location. 


The first three of these postulated sources of error were 
readily discarded. The calibration signals recorded with 
the strain gage signals were precisely those set into the 
Peeeeem prior to each shot; thus ruling out the possibility 
that either the calibration levelS were incorrect or that 
the playback amplifiers were out of calibration. To check 
the second point, a series of calibrated sine waves were 
input to the Fourier Analyzer using the same connections and 
equipment settings used for the strain gage data recordings. 
These calibrated signals were recorded on disk and then 
Sutput to the plotter in exactly the same manner as the 
Strain gage signals had been. In no case did the Fourier 
Analyzer alter the calibration of the signal. 

Whether the remaining suggested cause was, indeed, the 
source of error or whether the actual cause iS something 
Pamerely different cannot be determined at this point. 
Whatever the actual cause, the effect remains the same; the 
Strain magnitudes derived from the two underwater tests 
cannot be used for comparison with the results predicted by 
the USA-STAGS code. 

Although the magnitudes of the test derived strain data 
are clearly unuseable, the general shape of the strain 
curves should be essentially unaffected. iimecahiererempt CO 


/meamrrect the comparative nature of this study, it was 


ear 


determined that the shape of the curves from the tests 
Should be contrasted with the strain curves predicted 0 me 
USA-STAGS code. This attempt, too, was doomed to 

fLUStrar dome 

Figures 6.4, 6.6, and 6.8 are representative samples of 
the strain data calculated by USA-S7?4G5s functionine ae 
post-shot analyzer of the tee-stiffened plate. Figures 6.3, 
6.7, and 6.9 are the ceorrespondins St asec ce outputs.” 

Similarly, Figures 6.10, 6.12, and 6.14 are USA=5i3G. 
outputs for the rectangular-stiffened plate; and Figures 
6.11, 6.13, and 6.15 are their correspondine eri ain 
histories. 

While Similarities exist at some points between the 
actual strain histories and the USA~-STAGS predicted results, 
in general there is oe little poe ees This -was 
unexpected since earl ier tee c neers With personnel at 
Lockheed had indicated that extremely high degrees of corre- 
lation had been obtained when the code results had been 
compared to with tests conducted by the Defense Nuclear 
Agency and others. 


*2 unfortunately, after all of the problems encountered 


with the strain gages, there were not nearly as many of 
these aS might have been desired. 

aan again it should be emphasized that the shapes 
alone of the curves can be compared. AS already explained, 
the magnitudes of the strain gage records are not correct. 
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Figure 6.4 Tee-Stiffened Plate, Gage 2 (USA-STAGS ) 
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Figure 6.5 Tee-Stiffened Plate, Gage 2 
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Figure 6.8 Tee-Stiffened Plate, Gage 11 (USA-STAGS) 
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Figure 6.9 Tee-Stiffened Plate, Gage lil 
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Figure 6.10 Rectangular-Stiffened Plate Gage 1 (USA-STAGS) 
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Figure 6.11 Rectangular-Stiffened Plate, Gage 1 
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Figure 6.12 Rectangular-Stiffened Plate, Gage 15 (USA-STAGS ) 
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Figure 6.13 Rectangular-Stiffened Plate, Gage 15 
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Figure 6.14 Rectangular-Stiffened Plate, Gage 18 (USA-STAGS) 
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Figure 6.195 Rectangular-Stiffened Plate, Gage 18 
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we Unt oOmmeoss OMe Soluyions exist for this perform- 
pee, Of COUrse; including the extremes of either USA-STAGS 
Beeroe Strain histories being totally wrong. In this case, 
however, utilizing the extended capabilities of the code 
itself may have provided the most plausible solution. 

Mie USA postprocessor , POSTPR, has bie rea oa) iia OT 
m@epeculactcine the fluid pressures at the fluid structural 
interface. These results may then be displayed in either 
mpoular or graphical ee When this was done for the 
cases studied here, a very interesting result was obtained. 

Figures 6.16 and 6.17 represent the predicted fluid 
pressures at the center of the tee-stiffened and rectangular- 
stiffened plated, respectively. As can be seen, USA-STAGS 
predicts that intense negative pressures will be generated 
feeraid a2 very short time after the arrival of the shock 
wave. That is, USA-STAGS is predicting ee oe lO @au 
cavitation across the faces of these eon 

Viewed in that light, the reason for the seemingly OO ois 
performance of USA-STAGS in predicting the strain responses 


as noted earlier, no graphical capability exists on 


the VAX/VMS system at the Naval Postgraduate School. These 
post-shot analysis problems were run at Lockheed Research, 


mero Alto. It was they who most generously provided the 
meconopanying plots. : 
a2 


iment p  rocal Gay itatton would be marked by near 
zero pressures as the water flashed into steam to fill the 
void created by the sudden displacement of the plate. Since 
USA-STAGS assumes a homogeneous medium this must necessarily 
be represented by negative pressures. 
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can be better understood, and even the strain gage responses 
themselves become more understandable. 

For example, in Figure 6.11 an initial compressive Straim 
can be observed as the plate reacts to the incident pres- 
sure. This does not occur immediately but is delayed some- 
what as the inertia of the plate is overcome. In this . 
instance, the initial compressive strain 1S mimicked almose 
exactly by the response predicted by USA-STAGS. The strains 
then return to zero or become tensile aS cavitation occurs. 
As would be expected, the responses predicted by the code 
are much more violent then those seen in the strain gage 
histories. This is so because they are activated by large 
values of assumed pressure rather than the actual zero pres- 
Sure that. woulda) occur #Wil nea wera ten 

The other plots seem to confirm this basic trend: 
Initial zero or compressive strain depending upon the loca- 
tion of the element/strain gage followed by a positive trend 
in the Straim nistory as Cayvatat lon develoveme 

The actual mechanisms involved in the results that seem 
to be indicated here are not fully understood. The implica- 
tion, however, is that the velocity of the plate as it is 
being deformed is much greater than that of the fluid. aA 


Jaan is suggested that Figures 6.9 and 6.15 also adhere 


to this basic pattern, bibewate SUbNee 2 6O pe) eee wee 
problem mentioned earlier. It is believed that these plots 
Should actually be inverted. 
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discontinuity is therefore being created in the fluid/ 


mmemetvural interiace. This diScOntTinulty is manifes 


ct 


ed by 
the pressure drop which causes the local cavitation. 

iijewsotrain Gaee NIistorles Seem te imoly that local cavi- 
memeron 1S indeed occurring. Whether that is in fact the 
@eoe, Or whether the cavitation predicted by USA-STAGS is 
Simply a function of the way that the plates are modelled is 
something that must still be determined, however. 

In these studies, the plate alone has been modelled. 
Mmiewstructure that is ''kKnown’ to the computer code can 


Mime toro De Conceived aS infinite in extent and must neces- 


Sarily bear the full brunt of the pressure wave. ieee tu 
memeey, Chis 1S not so. The test plate is finite in size and 
pressure relief will occur around its boundaries. Since the 


Mimemath Of the shock front is on the order of five feet for 
charges of the weight used in these tests (this is a func- 
tion of the speed of sound in the fluid and the pressure 
decay, it will vary with charge weight), the entire assembly 
will be enveloped by the pressure front. This means that, 
almost instantaneously after the incident pressure is felt 
by the test plate, the pressure will begin to relieve around 
the edge of the plate and a negative pressure will be 
experienced on the backs of the flanges, and, soon after, on 
the rear of the backing Structure.” Therefore the entire 

> SC opal interference patterns can also be expected. 


ios 


force of the incident pressure Wave may not, 10 \teee oe 
impinging upon the test plate; since, while that wave is 
"oushing'' on the front side of the Structure, tere ae 
“pulling” on Ene Veae 

Whether this decrease in force would be sufficient to 
eliminate cavitation 1S an open question. More testing must 
be done. For the USA-STAGS code to be able to oreqi ae thas 
however, it is not enough to model only the test panel. The 
entire assembly of test panel, milled edge, and backing 
structure must be included. If localveavitation 12). eee 
predicted, then additional action will be required because 
Doubly ASymptotic Approximation which 1S the basi seeeen 
the USA code is not valid where discontinuities such as 
Cavitation exist between the fluid and the sti uetu. 


surface’> [meet 


C. CONCLUSIONS AND SUGGESTIONS FORP PUCiRE SsrUre 

Clearly, problems encountered both in conducting undem= 
water testing and in using the USA-STAGS code have precluded 
any firm conclusions from being made about the predictive 


capabilities of the code. Thus, the primary objeet1ve vem 





this study has not been met. It is in the fulfillment of 


> "Toeinieed has developed a processor, USA-STAGS CFA 


(Cavitated Fluid Analyzer) which can be used to model the 
cavitated region. The process that must be followed is to 
model layers of fluid at increasine Interval=e from tie pie 
until no further cavitation has been = found st tien o> ae 
the standard DAA approximation may be applied. 


Oe 


the second objective presented in the introduction, that of 
Mmeevlaine culdelines for future study, that definite contri- 
Dutions can be made. 

In the realm of underwater explosion testing, some 
soecific problem areas were noted in the preceding chapters. 
Solutions to some of these (e.g., strain gage choice and 
attachment) were presented in the text of those chapters. 
Others of these problems are of a more fundamental nature 
and must be discussed separately. 

Pest eaesauiSiactory Dian Of testing cannot be devel- 
oped without an adequate source of appropriately sized 
explosive charges. The inability to match charge weight to 
test platform dimensions makes it impossible to use any 
computer code in its preshot capacity and impacts negatively 
Meom the reproduceability of test results. Either the 
charge welght must be standardized or the test platform 
dimensions must be altered to accommodate changes in charge 
weight. In either case, the most apparent solution seems to 
Mena Carefully orchestrated test plan with long lead times 
Memgere Cach test shot. If amatch is not possible, it is 
Suggested that the shot should probably be aborted. That 
foe Ole I1ntOrmation Was obtained from such mismatched 
@emeoitionsS in this study 1S considered to be more a charac- 
teristic of a program in its developmental stages than of 
well-conceived experimental practice. Until more experience 
Mmemeatied, lb Appears that the rationale dictating use of 


aor 


the eight pound TNT charge as the standard |[Ret- 3) is Stem 
valid. It is therefore recommended that that weight charge 
be Maintained 1f lew eoss see 

second, the lack of adequate strain gage data has been 
an incapacitating flaw in the test results. Whether the 
poor quality of that data has been the result of ey Ca 
the equipment set-up at the West Coast Shock Test Facility, 
the interface between that equipment and the strain gages on 
the test plate, or some other cause 2s .uneer ra a) ee 
certain is that the present arrangement of sending the test 
plate and backing structure to WCSF a few days before the 
test shot for attachment and then allowing only a few hours 
immediately before the shot for equipment testing is 
insufficient. Development of an in-house instrumentation 
capability at the Naval Postgraduate Soneoll is considered to 
be essential. This would permit the gages and electronic 
equipment to be integrated into a single system. Thorough 
testing could be cGonducted before the actual date oil en 
test; and, equally important, responsibility for the onteews 
of the test would teside i 2 Sine t eee fants tee 

Changes in the electronic instrumentation have already 
been recommended [Ref. 3]. Enactment of those changes would 
undoubtedly improve the data responses. A few other modifi- 
cations suggested by the results of this study are consid- 


ered equally imvoortant, however. 
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Powe Mono Cali bration esesistors, the strain 
gage bridge should be balanced using strain gages 
identical to those on the test plate. These should 

be attached to an aluminum plate in the terminal box. 
This has two effects: First, the balancing gages 
would be subjected to the same changes in temperature 
as the plate gages. Second, the balancing leg of 

the bridge would be of the same length as the measure- 
ment leg. Combined, these would make the bridge 

far less sensitive to environmental changes.’ aA 

zero balance could be maintained and calibration 
fluctuations would be less likely. 

A zero-time triggering signal should be Superimposed 
[MOOI fast rument. channels. This would permit an 
exact assessment of the time of initial plate response 
to the incident wave. 

Further shielding and grounding of the pressure gage 
cables should be attempted to prevent the saturation of 
the pressure gage amplifiers through currents induced 


De enawac. devtOonat Lon . 


From the results that have been explored in the 


feeeree INS Section ©f this Chapter it is clear that certain 


Mulestions need to be raised about the validity of the flat 


plate model that has been used here. These same questions 


Should be asked about the ability of USA-STAGS or any other 


©eode to handle that model under the conditions imposed. 
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In particular the question of whether the local Gavita= 
tion predicted by USA-SIAGS Dseacw ee occurring er is 
Simply a result of the choice of model must Bewade-- - 

It is recommended that this study should be redone, but 
with the backing structure included aS a part of The yoda. 
model. POSTPR could then be run to determine whether cavi- 
tation was still cCC@currings aeress etn] pea If cavitat fem 
was still being indicated, then either a new model would 
have to be found or the cavitated region would have to be 
modelled using 4a code Capable wor saand aia 

Whatever the findings of continued tests, it may have to 
be recognized that increasing the scope of the models 
Studied may put USA-STAGS beyond the reach of masters-level 
research at the Naval Postgraduate School. USA-STAGS is a 
diverse, highly capable code which has the reputations. 
high levels of accuracy. The price that 1S pald i0r suena 
capability, however, iS size and complexily. §Givene1 ae 
limited time and computer resources available to most 
Students at NPS, USA-STAGS may be beyond their capacity. Lt 
Such research is undertaken, it 1S most emphatically recom- 
mended that a large block of time, weeks, perhaps, be spent 
at Lockheed Research in Palo Alto, California learning the 


code and perfecting models. 
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APPENDIX A 


SOURCE bs Ts 


Gireaeoenotix CONLAInS the information listings supplied 
with the STAGS and USA source tapes. Section A is a full 
listing of the contents of the STAGSS4 tape along with 
mies ructions for mounting the tape, command files for 
ereating and linking executables, and sample command files 
[imeerunnane the code. Section B contains similar listings 
for the USA source tape excepting that no command files are 
listed. USA places the link commands on the tape itself for 
immediate access. Command files for running the modules are 
Mereieical in construction to those for STAGS and can be 


adapted from those in Section A. 


las 


THE STAGSS4 SOURGES 2 ea 
Germs lilt) BUGUSe elaae 
STASSC1 PROGRAM TARE 
 VMAX/VUMS OPERATING Srsiem 
1 TAPE FORMAT — 
9 ea Gr 
1600 BP 
Labeled (STAGS4) 
< TAPE CONTENTS 
File Name Type Dbeseript ion 
1 INFO. VAX Tie Table Gr cones .¢, procedures 
2 STAGSU FOR Source FORTRAM program3 for ST4CS utility leepoeng 
S BLICMF MAR Source MACRO program *cr STAGS ur1licy Lior 
4 JOBSTA. MAR Source MaCPO orogran tor STAGS utility Library. 
5 STAGSi FOR Source FORTRAN progtsms ForSTAGSL lLiosrary 
6 MOUMNTECe ren Source FURTRAN programs for the MOUNT elemen >: 
& STAGSS FUR Source FORTRAN prograns for STAGS2 Library 
8 Steril On Source FORTRAN prog-ams for ST4PL Library 
7 POST FOR Source . PORLRAN prog-am>. for. POS ie lvoe 
40 ST4C6S1,0L8° tLrersre Libraeyu or S256]. Fo coc can 
a STAGSe OLB Library C1bTrery of SiS272250eo7rams. 
12 STAGSU CES) Serb sare Library or S° “Ve OGriity rome. 
te POST. Ces Library Ligrary of PIelr prea ranis 
14 STAPLE. OEs Library Library of Site plo: onedn omer 
iS Nice OL LA bra gy 
16 Te AAPOR Source & ThUTs conen fo. occ on 
17 TPeOG Ob sect THUES COM eC roces>s cour. 
18 ZS7S. FOR Source Z-system utilities (for TP). 
LY RES. FOR Source Reduced Equation Solver 
20 RES OBJ Obyect Reduced Equation Solver 
at 235s Os Library Z-system utiizties object Library 
22-36 CASE#. *. # Last Input and ovtcut for so somo le eoo--. 
J; LOADING FROM TAPE 
Use the following control cards to (copy “all 71 less ceca 





$ MOUNT MTAO: STAGS4 TAPE 
$ COPY TAPE: *.# #* (Jj 


Files may also be read indidually 
following control cards ‘(e=gre 6 laac 


rTrom cape with the 
BlAGs wo eae 


$ COPY TAPE Sracsl OER £3 


Complete sample case output may be odtsired as roliows 


* PRINT HEAD CAse ss. = 
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em Ww AR t4 fa 


www th tH 


w 


=_—— -) 


sy Eten 2 Or WEOmciie oF os r ams ame ler 
Bees lAGo! FOUR7MILIS?T 
Berea Y/UOREA4TE SiaGsy OLS SsiAGsu 
MACRO BLIOMF MAR 
MACRQ JOBSTA MAR 
Lid STAGaY ELIOMF, UGb=TA 


ByeSToeS1 -OR/HOLI< 
LIBPARY/CREATE STAcs 
FOR MOUNTEL FOR/NOLI 
LiB STAGS1 MOUNTEL 


Bes SiaGs! 


fi 
¢ 
4 
Si 


gee siAGSe FOR/NOLIST 

MrgRearny/CREATE STACSS 

momemovel=t, FORSNOLIST 
ie) S'TAGSe MOUNTEL 


ee 


mere STAPL FOR/NOULISI 
MeeprersY/CREATeE STarPl SLB STAPL 


FOR FPOSTP FOR/NOLIST 
LIBRARY/CREATE POSTP OLE POSTP 
FOR TP FOR/NOLIST 

FOR ZSYS FOR/NOLIST 

LIB/CREATE ZSYS ZSYS 


gem mes. -GR/NOLIST 


B. Bei AGsi CoM (Load StAysl) 
LINK/MAPS=STAGSI/BRF Tere cf ECUTABLERHSTACSI- 
SraGol @es/eis7 INCLUDE=H4EIN1L, — 
STACSU. GLS/EIB 
on LSTAGS2. COM (Load STAGS2) 
LINK/MAP=STAGS2/ BRIEF /EKECUTABLE=STAGS2- 
STAaACSe. OLB/LIB/ INCLUDE=MAING, - 
STAGSU OLB/LIB 
D. Star eo (Load STAPL) 
LINK/MAP#STAPL/BRIEF/EXECUTABLE=STAPL- 
saree. Ceeyv el B/ INCLUDE =STarl, — 
SracsU, ULB7EL15: - 
(File containing 


lecal olot routines) 


es Pee See (Load FPUSTP? 

mete MAP=POSTP/ERIEF EXECUTABLE =POSTP- 
Bootes Oely ley hNECUDe = 5 ir. — 
Siebol OCE/L1e,SieCS2 CLB/Lib.— 
Brace GEE/LiesittlCe, OLB VierE& 
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Pee HR FW Wm HH 4 Rw wm ew ott Ub 


FARA HH 


E Sia so bec 

=SitN Cé&sei ie 
WSS isis ASS ewe 
S556 C4521 MOD 
RUN STAGS1 EXE 
PR tit © GASe Sua 

G. STAGS. COM 
ASSIGN CaAsei INP 
ASSIGN Castel. QUT 
ASsISWl CASSEL MOD 
aUN STAGS! EXE 
45S5IGMm CASE1 RST 
AZSIGN CASE! SOD 
RUN STAGSS EXE 


PRINT CAsei gui 


H.. -Pasie Teer 


ASS boN@e aS = 4c Lie 
ASSIGN CASES OUT 
RUM Sega eke 


ASSIGN 
ASSIGN 


CA5E4 PP 
CASc- FUP 
ASSIGN CASc4 RST 
ASSIGN CAsc4 SOD 
A oe oie 4 
BEL -F GR@GOs- DAT 


1 Si APG. COM 


ASSIGN CASet PIN 
ASSIGN G4sel RST 
ASSIGN CASEL MOD 
ASSIGN CASE. QUT 
RUN STAPL. EXE 


T 29s" acne 
2a ved" 


POR aS 
See 
FOrtst 

(21, e@cCut] St7AGsl tance =) 
SOeiwe = 
FORUGS6 
FORV21 (eptional if Mode. Clots os eo 
FORGEG (previous solwUtiaor gata, 3s" don 
rORGE2 (spbional if selucisn cesamseog 1c orc 

(Execute STAGSI and POS; 
FORGGS 
FOROG6 
(May bee omitted iy FORGQ] bam scaved: 

FHRUOS 
FORGGS 
FORD20 (previous solustian data from CASES SOD) 
FORGeE2 


(Execute STA&FPL) 


FOROQOS 
FORGE 

Penge) 
FOROOS 


(Execute plot file on local plat device) 
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B. THE USA SOURCE TAPE 


scene er Srl NG 
WAsA TAPE LOUACING INa 
Meee a0 Cmeeat ay As “Suma si1=s+ trPemue = 
Peeeuleaglos, and manuai:s -or 32 "umoer or Us-. 
meemecneatec with COPY commanes and can be - 
Tereoaiy a -ew Lines of input at tne fermi: 
directe-y will reveal tnat tne Fries nave -s 
Meomconme “ey to keeping tnose ftiies tagether - 
mpeeeiassct. ihe foliaqwinge table conteins this 
USe eo = nweaS=-e- 
Bernas 
AUGMAT 
Pris 
eS 
MODE se 
fomergag the @tIMINT files, for exampis, 30 =: 
preomero restde and simoly do 6 COPY TAFE # »* 
the incternai name that has oeen assigned 3 
Meeress 15 simoly repeated umti!l ell aesire: 
Loaded 
Moye Ss character erteneciom Fon tne =. 
mmeaceording £9 the fellewing table 
Sx, Tes! ON 
ee UP 
NG 
. VAX 
7OCE 
ecour 
Ree 
. VAL 
ror example, AUGMAN UAL is a File containing 
MnO SsAS. COM is the LINK statement for USA-s 
that all LINK statements must be modified s. 
referenced in them retlect names that are in 
Bom@ieaim wsers mames that will not exist at 9 
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ay 
ip “b +B 'D 
sad 


my 


- 


mv wa 


Woe. cord Wie 
Srnand whera 
7ape 


Rave 


Hayy anes 


INCLUDE PROCEDURES 
INCLUDE PROCEDURES 
ES TrMeol Cs 

OBJECT LIGRARIES 
LINK STATEMENTS 
EXECUTABLES 
MANUALS 


eEhe AUGMAT manual whrile 


TAGS feeenouwld oe nosed 
1olecrUeeais che cirectories 
use here at Lockheed and 
ther VAs installations 


Dich Ss Geeee eee 
Almass Sxe. t 37/57 e-SES-195. 07 0 
ly ee | Ae Deed 17 LeaeMaAr HL 95— 1 OD 
AMSEO! FUR. i ed o= JUNNW TSB. 1) OO 
aloe0a FOE im = se = NIN a oe eee) 
AMGEO3 FOP. 1 mS i 19--JGE = Ss 2 oe ga 
ANGEO4 PDP, 1 7 ae i= JUN= 1 SSeS eeenow 
PJP ese n AG/A9 S-SEF-1984 oC 30 
FLUMAN UAL: 1 a ae S—SEP-f984 35° VC 
PLUMAS FOR: tl 74/74 a SePHoe  eee U 
LIMNAFLU, COM. 4 ee E—SEP-198= 7 GO 
NGMFELU PDP: 1 17 Pe -APR~-1SBL 2 GO 
AMGEm tm terse 7} 4—-JUN-155- 00 GO 
AMGeOS FDPae leash eo -JUL = 1 Sea eo 
ACES Oe 25/25 0° See e231 
AUSMAN UA e Ley a B—MAY—v¥se0 we. Co 
AVWGMAT FOR, 2 acy &e-Ser-i:96- 109 Go 
FILOGAT POP: 2 ty 23-Mar-i S84 Ue GO 
| INKAUG COM, 2 | é—-SEP-1i9e8- vi. OG 
NAME G ee ec ase - eyed SenmOPR-iSS2 29 GO 
PREFROC Exe. 2 wee &-SEF-i198- 70.60 
TRAWSE FOP: 2 1a 2F-MNAY~-i 98s 30° OG 
GrHisean ee Ona Ly 30-MAR-i93- Ge: oO 
COMaDE, FE De: 2 7 BO-MAR-1L95= oo OO 
JSKFAC. INC, 3 ee 2O—-MAR-—178e UG OO 
DATE Li eee, s ey} 15-JUN-1964 93 GO 
DATINT PDP; 3 tee Wee BO-MAR-1 98s Yt OO 
LATTE Oleree S27 Sa 6-SEP-198+4 02:00 
INTMAN VAL, 3 ea/e2 3-NOV-1982 20 60 
LINKUSA COM. 3 P71 &-SEP-198+ O80 OO 
N&MINT. PDP; 3 7 SO-MAR-1984 ul OO 
TIMINT FOR: 3 $8758 6-SEP~-1984 OC OO 
UNWASH EXE: 3 70/70 &-SEP-1984 00: 00 
COMADD. PDF: 4 fad SO-M4R-198s GO. OO 
DATFIL. PDP; 4 Ly 3O-MAR-1984 C60: 00 
LINKPOS. COM: 4 1/4 3O0-MAR-1984 00: 00 
NAMINT. PDP: 4 Lv S3O-MAR-1984 GO: 00 
POSLIB. OLB. 4 36/38 30-MAR-19843 G0 6O 
POSMAON UAL: 4 Loft e7-MAR-1998= OC. OO 
POSF-RO FOR; 4 41/41 OG-MAR-iFS- UC. 00 
POSTPR. EXE, 4 61/61 B3O—-MAR—-1984 Go: GO 
CHGDAT. PDP; $ 17a S-APR-1963 00 OO 
COMADD. PDP; § 7 18-MAR-19S3 00: 00 
DATFIL. PDP, S Pe 15-JUN-1984 O00. 00 
D&TINT PDP: S Ld 15-JUL-1982 OC GO 
LINKUSAS. COM; 5 : Lt 64-SEP-1985 00 GOO 
MODLIS Pek Loy ae 6-SEP-19584 00:00 
MGCL IS Wie, G/9 &—SEP~198< OC OO 
NGQMINT. PDP: > Let ea-APR-1982 G0O° 00 
UeAc Sten ae . eeye22 4-SEP-1981 99 00 
UST. 12 oc. 5  * Se Vi 45-SEP-1984 GC GO 
US TSrM er Ons > Lee O~SEP~-i98s OC CO 
RBiOQSUP MAR. 8 - O79 15-4UG-198i “CG GO 
CC FOR; 8 171 S-SeP-198l 9 50 
GMERR. PRC. 3 lee H-JAN-19S5S Jo CO 
CitsaSP FOR, 8 La/ig BG-JAN- 1982S 9G GO 
DMGAaSP OBU: 5 eae G-GAm-19SBl OC O00 
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APPENDIX B 


BASIC SHELL TYPES Cro) ith >t aS 


This appendix contains a detailed listing of ~ he wele veg 
pane Shell types used with the STAGS code. Drawings are 
included to show numbering conventions for edges and 
corners. This information Gam be particularly Weer wee 
discretizing the shell unit and specifying compati>oi 


with ,adjoining shells. (All figures from Reterence mia 
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USa=STAGS INPUT ares 


This appendix contains the DQap@ert 1 lessee ecu] 
complete USA-STAGS runs for both the tee-stiffened and 
rectangular-stiffened plates. These are listed in the order 
in which they would normally be run: 

- STAGS2 (Modal Analysis) Input File 

_ STAGS1 inoue soe 

- Fluid Mass Input 

- Augmented Matrix Input 

- Time Integrator Input 

- USA postprocessor (POSTPR) Input 
The input files for the tee-stiffened plate include an 
example of STAPL input. The files for the rectangular-— 
stiffened plate contain an example of PUSteei ip Tice 
Since these last files are readily generalized, it was not 


considered necessary to provide examples for each case. 
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